Absorption of mass and angular momentum by a black hole: Time-domain formalisms 
for gravitational perturbations, and the small-hole/slow-motion approximation 
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The first objective of this work is to obtain practical prescriptions to calculate the absorption 
of mass and angular momentum by a black hole when external processes produce gravitational 
radiation. These prescriptions are formulated in the time domain (in contrast with the frequency- 
domain formalism of Teukolsky and Press) within the framework of black-hole perturbation theory. 
Two such prescriptions are presented. The first is based on the Teukolsky equation and it applies to 
general (rotating) black holes. The second is based on the Regge- Wheeler and Zerilli equations and 
it applies to nonrotating black holes. The second objective of this work is to apply the time-domain 
absorption formalisms to situations in which the black hole is either small or slowly moving; the mass 
of the black hole is then assumed to be much smaller than the radius of curvature of the external 
spacetime in which the hole moves. In the context of this small-hole/slow- motion approximation, 
the equations of black-hole perturbation theory can be solved analytically, and explicit expressions 
can be obtained for the absorption of mass and angular momentum. The changes in the black- 
hole parameters can then be understood in terms of an interaction between the tidal gravitational 
fields supplied by the external universe and the hole's tidally-induced mass and current quadrupole 
moments. For a nonrotating black hole the quadrupole moments are proportional to the rate of 
change of the tidal fields on the hole's world line. For a rotating black hole they are proportional 
to the tidal fields themselves. When placed in identical environments, a rotating black hole absorbs 
more energy and angular momentum than a nonrotating black hole. 



I. INTRODUCTION AND SUMMARY 
A. Goals and motivations 

The work described in this article is concerned with 
the absorption of energy and angular momentum by a 
black hole when physical processes in its exterior produce 
gravitational radiation. It is assumed throughout that 
the rates of change of mass and angular momentum are 
sufficiently low that they can be calculated within the 
framework of first-order perturbation theory, in which 
the black hole differs only slightly from a stationary and 
axisymmetric Kerr hole. 

The first goal of this work is to obtain practical pre- 
scriptions to calculate the black-hole absorption, and to 
modernize the tools fashioned in the early seventies by 
Teukolsky and Press lj . An essential aspect of the new 
prescriptions is that they present the absorption formulae 
in the time domain instead of the frequency domain; they 
presuppose that in accordance with current trends, the 
equations of black-hole perturbation theory have been 
solved as partial differential equations in the time domain 
instead of ordinary differential equations in the frequency 
domain. Two such prescriptions are presented here: the 
first is based on the Teukolsky equation 2] and it ap- 
plies to general (rotating) black holes, while the second 
is based on the Regge- Wheeler 3] and Zerilli 4] equa- 
tions and it applies to nonrotating black holes. 

That the Teukolsky equation can be separated in all of 
its variables is surely one of its most important proper- 
ties. To a large extent, it is this property that has permit- 
ted progress during the continuing exploration of physical 
processes taking place in black-hole spacetimes (see, for 



example, the book by Frolov and Novikov, Ref. [5|]). But 
it has to be acknowledged that the historical importance 
of the separation property has diminished in recent years, 
as a number of time-domain integrators of the Teukolsky 
equation have been developed [y, H> & an d P u t to 
use in various applications. The numerical task of solving 
the Teukolsky equation in the time domain is still chal- 
lenging: after decomposition into azimuthal modes one 
must solve for a function of time and two spatial coor- 
dinates. But time-domain methods appear now to be at 
least competitive with frequency-domain methods, with 
which one must solve for a number of radial and angular 
functions, the number increasing as the spectrum of rel- 
evant frequencies becomes wider. And it appears likely 
that the future will witness an increasing dominance of 
time-domain methods over frequency-domain methods. 

While the superiority of time-domain methods is still 
to be proved in the case of the Teuko lsky equation, it has 
clearly been established [ITIJ12I IT3L Il4l ITal in the con- 
text of the Regge- Wheeler y| and Zerilli [4J equations, 
which determine the metric perturbations of a nonrotat- 
ing black hole. In these cases the angular dependence of 
the perturbation variables can be completely separated, 
and the integrator faces the relatively simple task of solv- 
ing for a function of two variables (time and a radial 
coordinate). Simple, but powerful, numerical methods 
have been devised for such problems, and these can 
even handle, without approximations, a singular source 
term contributed by a point particle. The time-domain 
Regge- Wheeler and Zerilli equations are thus very easy 
to integrate, and there is now little reason to go back to 
a frequency-domain formulation. 

The new popularity of time-domain methods to solve 
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the equations of black-hole perturbation theory calls for 
new prescriptions to calculate the black-hole absorption 
of energy and angular momentum. The only recipe cur- 
rently available is the formalism of Teukolsky and Press 
Q , which is based on the frequency-domain formulation 
of the Teukolsky equation Q • This formalism is not well 
adapted to time-domain calculations, and in this work I 
provide the required translation of the Teukolsky-Press 
recipe to the time domain. Another limitation of the 
Teukolsky-Press formalism is that although it can be 
applied without difficulty to a nonrotating black hole, 
this requires the use of the Teukolsky equation instead 
of the more practical Regge- Wheeler and Zerilli equa- 
tions. Another objective of this work is therefore to re- 
late the absorption of mass and angular momentum by 
a Schwarzschild black hole to the time-domain solutions 
to these equations. 

In effect, this work is about providing practical time- 
domain formulae for the fluxes of mass and angular mo- 
mentum across a perturbed black-hole horizon. For a 
nonrotating black hole these formulae are based on the 
Regge- Wheeler and Zerilli equations, which govern the 
behavior of the metric perturbations. For a rotating 
black hole the formulae are based instead on the Teukol- 
sky equation, which determines the perturbations of the 
Weyl curvature tensor. 

The second goal of this work is to apply the time- 
domain absorption formalisms to physical situations in 
which the black hole can be considered to be either small 
or slowly moving. In the context of this small- hole/slow- 
motion approximation (which I will describe in Sec. I E 
below), the equations of black-hole perturbation theory 
can be solved analytically, and explicit expressions can be 
obtained for the absorption of mass and angular momen- 
tum. While many results have been obtained along those 
lines in the past ll 71 Il8l Il9| , they were all restricted 
to various special cases; the results presented here con- 
solidate and generalize these previous works. 

The absorption of mass and angular momentum by a 
black hole is generally very small. In particular, the effect 
is likely to be too small to be observed in a gravitational- 
wave signal that would be measured by ground-based 
detectors such as LIGO, VIRGO, and GEO600. For ex- 
ample, Alvi 01 has calculated that for binary systems 
involving holes with masses ranging from 5 to 50 so- 
lar masses, black-hole absorption is truly negligible: It 
contributes only a small fraction of a wave cycle during 
the signal's sweep through the detector's frequency band. 
For this type of source the tools developed in this paper 
are not needed. 

In some circumstances, however, the black- hole ab- 
sorption is a significant effect that should not be ne- 
glected . In particular, it is likely to be observed 
in gravitational-wave signals that would be measured by 
a space-based detector such as LISA. For example, Mar- 
tel ^3 has shown that during a close encounter between 
a massive black hole and a compact body (of a much 
smaller mass), up to approximately five percent of the 



total radiated energy is absorbed by the black hole, the 
rest being transported out to infinity. Hughes [2l| has 
calculated that when the massive hole is rapidly rotating, 
the absorption has the effect of slowing down the inspiral 
of the orbiting body, thereby increasing the duration of 
the gravitational-wave signal. For example, a 1 Mq com- 
pact body on a slightly inclined, circular orbit around a 
10 6 Mq black hole of near-maximum spin would spend 
approximately two years in the LISA frequency band be- 
fore its final plunge into the hole; Hughes shows that the 
black-hole absorption contributes approximately 20 days 
(and 10 4 wave cycles) to these two years. For this kind 
of situation the absorption is important, and the tools 
developed in this paper will be useful. 



B. Perturbative methods 

A natural starting point for the calculation of black- 
hole absorption would be the definition of a dynamical 
mass M(v) and angular momentum J(v) on a cross sec- 
tion v — constant of an evolving event horizon; here 
v is a suitable advanced-time coordinate on the hori- 
zon. Armed with such definitions, one would differen- 
tiate with respect to v to obtain M(v) and J(v), and 
seek to express the right-hand sides in terms of standard 
perturbation variables. Such an approach to black- hole 
absorption has recently been pursued by a number of 
workers EHHIHIHl, and the resulting (inequivalent) 
formalisms can be formulated exactly in fully nonlinear 
general relativity. These formalisms are based not on the 
event horizon, but instead on the hole's trapping horizon, 
a generally spacelike hypersurface foliated by marginally 
trapped surfaces; and in the Ashtekar-Krishnan formal- 
ism I22II23I] the definitions for M(v) and J(v) come from 
the Hamiltonian formulation of general relativity. These 
formalisms are interesting (and useful in the context of 
numerical relativity) because they are fully general, and 
because they involve a hypersurface (the trapping hori- 
zon) whose intersection with a given Cauchy slice is easy 
to identify; the event horizon, on the other hand, can be 
identified only once the future history of the spacetime 
is completely known. 

The approach adopted here to calculate the black- 
hole absorption is not the one described in the preceding 
paragraph; it is based instead on black-hole perturba- 
tion theory, and it assumes that the evolving black hole 
is only slightly different from a stationary and axisym- 
metric Kerr hole. Because the analysis is restricted to 
first-order perturbation theory, it is possible to proceed 
without the specification of a mass function M(v) and 
an angular-momentum function J(v), so long as only the 
long-term changes in mass and angular momentum need 
to be calculated. In this long-term view one imagines 
that the black hole starts in an initial stationary state 
characterized by the parameters (M, J) , is perturbed for 
a time Av by some external process, and then returns 
to another stationary state characterized by the param- 
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eters (M + SM,J + 5 J). One then defines the aver- 
aged rates of change of mass and angular momentum by 
(M) = (6M)/(Av) and (J) = (SJ)/(Av), and one ma- 
nipulates the equations of black- hole perturbation theory 
to calculate these quantities. This is what I set out to do 
in this work. The perturbative techniques demand that 
SM <C M, 5 J <§; J, and the long-term view demands 
that Av 3> M. While the price to pay is a substantial 
loss of generality with respect to an exact formulation, 
the perturbative-long-view approach adopted in this pa- 
per allows one to proceed without having to choose a 
specification of M(v) and J(v), with the derived benefit 
that the final results are robust with respect to a change 
of definitions. Another benefit is that the approach is 
based on the event horizon (the true boundary of the 
black-hole region) instead of the trapping horizon; while 
locating the event horizon in a nonlinear situation can be 
difficult, this is not a problem in the perturbative-long- 
view approach. 

The mathematical techniques required for the descrip- 
tion of a perturbed event horizon go b ack to the pio- 
neering work of Hawking and Hartlc 26] , and these form 
the basis of this work (they also formed the basis of the 
Teukolsky-Press prescription p|). These techniques are 
reviewed in Sees. II, III, and IV of the paper. I begin in 
Sec. II with a description of the unperturbed horizon of 
a stationary Kerr black hole. In Sec. Ill, I consider the 
dynamics of a general evolving horizon, and in Sec. IV, I 
specialize the discussion to event horizons that are per- 
turbed versions of the Kerr horizon. The equations that 
govern the behavior of the horizon's null generators are 
given in Sec. Ill and IV in a form that closely resembles 
the treatment provided by Price and Thorne [2?|. and 
in Chapter VI of the book by Thorne, Price, and Mac- 
donald |2^|. Although these equations are well known, 
I derive them ah initio in order to clearly identify the 
simplifying assumptions that are incorporated along the 
way; in particular, the averaging procedure involved in 
calculating (M) and (J) is explained fully in Sees. IV B 
and C. 

The main ideas behind the Hawking-Hartle techniques 
|26j are as follows. The intrinsic geometry of an evolving 
event horizon is described by a two-dimensional metric 
tensor jab, which depends on the advanced-time coordi- 
nate v as well as angular coordinates A (A — 2, 3). The 
metric is degenerate (and explicitly two-dimensional) be- 
cause the horizon is a null hypersurface; v is a parameter 
on the horizon's null generators, and A are generator 
labels that stay constant as the generators move. The 
evolution of jab is determined by the behavior of the 
generators, which is described in terms of an expansion 
scalar 8 and a shear tensor (Tab- The evolution of the 
shear is driven by the spacetime's Weyl curvature, and 
the evolution of the expansion is driven by the square of 
the shear tensor. The evolution equations can all be in- 
tegrated (assuming that the Weyl curvature is specified), 
and the solution for (0) determines (A), the (averaged) 
rate of change of the black hole's surface area. The final 



step is to use the first law of black-hole mechanics (see, 
for example. Chapter 12 of Ref. 29], or Chapter 5 of 
Ref. [3(j ) to relate this to (M) and (j). The end result is 
Eqs. KT^ - KT^k in Sec. IV C, which express (M) , (J), 
and (A) in terms of o~ab contracted with derivatives of 
jab, integrated over a cross section of the event horizon. 
These equations are not new: they were first presented 
by Thorne, Price, and Macdonald |28|. but the derivation 
provided here is substantially different from theirs. And 
while these equations are not themselves very practical, 
they form an excellent starting point for the development 
of practical formalisms. 

C. Curvature formalism 

The development of a time-domain formalism to cal- 
culate (M), (J), and (A) in terms of standard curvature 
variables is undertaken in Sec. V. The end result of this 
reformulation of Eqs. H4.22fl - 14.24fl is the following pre- 
scription (Sec. V D) to calculate the black-hole absorp- 
tion. 

First, define a Teukolsky function W = — i/>o(HH) as in 
Eq. (|5.6fl , in terms of a null-tetrad decomposition of the 
perturbed Weyl tensor. The label "HH" indicates that 
the Weyl tensor is decomposed in the Hartle-Hawking 
null tetrad j2|J , which is well behaved on the future hori- 
zon of the Kerr spacetime. 

Second, decompose the Teukolsky function in terms of 
azimuthal modes proportional to e™^, 

oo 

*(v,r,8,^)= * m (v,r,8)e^, (1.1) 

m— — oo 

where m is an integer. Because the Kerr spacetime is 
axially symmetric, each mode ty m (v,r 1 0) evolves inde- 
pendently. Note that the coordinates (v,r,0,ip) are in- 
going Kerr coordinates (Sec. II A), and that they are well 
behaved on the event horizon. 

Third, integrate the Teukolsky equation Q for each 
relevant mode ^ m (v, r, 0), and evaluate the result at r = 
r + = M + V M 2 — a 2 , the position of the unperturbed 
horizon; a = J/M is the specific angular momentum of 
the Kerr black hole. 

Fourth, calculate the integrated curvatures 

/•OO 

$ r £(v ) 6) = e KV e - ( - K ~ mnH ' )v 'y m (v',r + ,8)dv' (1.2) 

J V 

and 

$™(v,0) = j e imQsV '^ m (v',r + ,0)dv', (1.3) 

J —oo 

where n = (r + — M)/(ri + a 2 ) is the surface gravity 
of the Kerr horizon, and £1h = + a 2 ) its angular 

velocity. Notice that $™ at advanced time v depends on 
the behavior of ^ m at later times] this is a consequence 
of the teleological nature of the event horizon. 
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Fifth, and finally, insert the integrated curvatures and 
their complex conjugates (indicated with an overbar) into 
the flux formulae 



(M) 



r\ + a 2 



4k 



E 



m— — oo 



2k / (|$™| 2 )sin(9d6> 



imO H / ($™$ m - $™$ m ) sin 9 d9 

(J) = --^ E (*») 



(1.4) 



4k 



771 — — OO 



and 



K 

8^- 



771 — — OO 



(|$™| 2 )sin6»d6». (1.6) 



These equations reduce to those of Teukolsky and Press 
when t> m {v,r,9) is a pure mode of frequency u>, 
^ m oc e -4 ""'; this is established in Sec. V C. Equations 
(|1.4(l - l|1.6|l are therefore the time-domain equivalent to 
the standard frequency-domain prescription. 



D. Metric formalism 

The curvature formalism of the preceding subsection 
applies to a general rotating black hole, and the special 
case of a nonrotating hole can be handled simply by set- 
ting a = 0. But in this case it is often desirable to work 
with metric perturbations instead of curvature perturba- 
tions, and it becomes useful to present the flux formu- 
lae in terms of ^ l ^{v,r) and ^^(w, r), the standard 
Regge- Wheeler and Zerilli-Moncrief functions, instead of 
the Teukolsky function *5> m (v, r, 9). Here the decomposi- 
tion into modes involves spherical-harmonic functions of 
degree / and azimuthal number m. 

The development of a time-domain formalism to cal- 
culate (M), (J), and (A) in terms of standard metric 
variables is undertaken in Sec. VII, after laying some im- 
portant foundations in Sec. VI. The end result of this 
reformulation of Eqs. ((4.22|) - ((4.24() is the following pre- 
scription (Sec. VII C) to calculate the absorption of mass 
and angular momentum by a Schwarzschild black hole. 

First, integrate the Regge- Wheeler equation Q for 
^Rw( w ' r )' which describes the odd-parity sector of the 
metric perturbations. This gauge-invariant function is 
defined in subsection 3 of the Appendix. 

Second, integrate the Zerilli equation for ^2m(v, r), 
which describes the even-parity sector of the metric per- 
turbations. This gauge-invariant function is defined in 
subsection 4 of the Appendix. 

Third, and finally, evaluate the Regge- Wheeler and 
Zerilli-Moncrief functions at r = r + = 2M and insert 



them into the flux formulae 

^E E (i-w+w+2) 



(M) = 



1=2 m=-l 



and 



(J) = 



x (4|*^(z;,r + )| 2 + |^( W ,r + )| 2 ) (1.7) 

^ OO I 

^EE (i-mi+m+2)(im) 

1=2 m = -l 

x (i^( V ,r+) J V ¥^(v',r + )dv' 

+ ¥iUv,r+)^(v,r + )). (1.8) 



Except for the substitution (v — > it, r + — > oo), these for- 
mulae are identical to Eqs. (|A.26|I and (|A.27|) which give 
the rates at which energy and angular momentum are 
transported to future null infinity. Note that for a non- 
rotating black hole, the first law of black-hole mechanics 
reduces to (k/8tt)(1) = (M). 

The flux formulae of Eqs. (11.7(1 . (|1.8|) were first pre- 
sented and used by Martel [Ij] in his numerical explo- 
ration of gravitational-wave processes associated with 
the motion of a small-mass body in the field of a 
Schwarzschild black hole. Although he arrived at the 
correct results, the derivation of Eqs. 1(1. 7J) and 1(1.8(1 pre- 
sented by Martel is flawed, and the analysis presented in 
Sec. VII puts them on a firm footing. Martel's derivation 
incorporates both a conceptual and a computational er- 
ror, the latter compensating for the former. Martel based 
his derivation of Eq. ((1.7(1 and 1(1.8(1 on Isaacson's effec- 
tive stress-energy tensor for gravitational waves [3ll l32| , 
incorrectly assuming that Isaacson's high-frequency de- 
scription is always applicable near the event horizon of a 
black hole. This story is related more fully in Sec. VII 
C, and its proper telling requires the connection between 
(M), (J) and the Isaacson stress-energy tensor estab- 
lished in Sec. VI C. The limitations of the high-frequency 
description become especially clear in view of this con- 
nection. 



E. Small- hole/slow-motion approximation 

A concrete evaluation of the flux formulae would typ- 
ically require the numerical integration of the Teukol- 
sky equation, or the Regge- Wheeler and Zerilli equations; 
an illustration is provided by Martel's recent work [l4| . 
But in some circumstances it is possible to solve these 
equations analytically, and to obtain approximate expres- 
sions for (M) and (J). I carry out such calculations in 
Sees. VIII and IX, in the context of a small- hole/slow- 
motion approximation that I now describe. 

Consider a situation in which the black hole is im- 
mersed in an external universe whose radius of curvature 
1Z is such that M/1Z <C 1. For example, suppose that the 
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black hole is moving on a circular orbit of radius b in the 
gravitational field of another body of mass M oxt . Then 
TlT 1 is of the order of the hole's angular velocity, and we 
have 



M M 

n ~ m + A/ ox t 



r 



M + M r 



where V is the hole's orbital velocity. One way to make 
this ratio small is to let M/M cxt < 1; then M/TZ will 
be small irrespective of the magnitude of V. This is the 
small-hole approximation, which allows the small black 
hole to move at relativistic speeds in the strong gravi- 
tational field of the external body. Another way is to 
let V -C 1; then M/TZ will be small for all mass ra- 
tios. This is the slow-motion approximation, which al- 
lows the slowly-moving black hole to have a mass com- 
parable to (or even much larger than) M cxt . These two 
limiting approximations are special cases of the funda- 
mental requirement that M/TZ be small; I call this the 
small- hole/slow- motion (SH/SM) approximation. 

When viewed on the large scale TZ, the black hole occu- 
pies a very small region of the actual spacctime, and this 
region can be idealized as a world line 7 in the external 
spacetime. Let u a be the (normalized) tangent vector to 
this world line, and call this the four-velocity of the black 
hole in the external spacetime. Assume that the Ricci 
tensor of the external spacetime vanishes on 7, so that 
the black hole's neighborhood will be empty of matter. 
The curvature of the external spacctime in this neigh- 
borhood is then described entirely by the Weyl tensor. 
This can be decomposed into its electric and magnetic 
components (see, for example, Ref. H3), respectively 



£a/3 = C lliav f}U tl U V 



and 



B a p = -u'V'C 



(1.9) 



(1.10) 



where the Levi-Civita tensor e mv p and the Weyl tensor 
Cfiaup are evaluated on the world line 7. The tensors £ a p 
and B a j3 are orthogonal to u a , and they are both sym- 
metric and tracefree; they comprise all ten independent 
components of the Weyl tensor. These tensors represent 
the tidal gravitational fields that are supplied by the ex- 
ternal universe, and these act on the black hole so as to 
produce a tidal distortion. This distortion, in turn, gives 
rise to a change of mass and angular momentum that 
can be computed with the formalisms described in the 
preceding subsections. 

In Sec. VIII, I calculate (M) and (J) for a 
Schwarzschild black hole moving in an external universe, 
to leading order in a SH/SM approximation. The results 
are 



(M) 



16M e 



45 



{Lp£ a f + B a pB afi ) 



(1.11) 



and 



(J) 



32M 6 
45 



u"e tienS (e a p £*> + B%B^)s 5 , (1.12) 



where s a is a unit vector, orthogonal to u a , that gives 
the direction of the vector (J a ) — (J)s a (a more precise 
definition is found in Sec. VIII F), and £ a p = £ a /3;txU >1 , 
B a p = Bafcuii 1 * are the proper-time derivative of the tidal 
gravitational fields. From these expressions we infer that 
(M) scales as M 6 /TZ 6 , while (j) scales as M 6 /TZ 5 . In 
Sec. VIII G, I show that the change in mass and an- 
gular momentum can be understood in terms of a cou- 
pling between the tidal fields and the hole's induced mass 
and current quadrupole moments, which are given by 
M aP = f§M 6 £ Q/3 and J afj = ^ M & B a p., re spectively. As 
illustrative examples, Eq. (|1.11|) and l|1.12[l are evaluated 
in two different limits in the case of circular binary mo- 
tion: In Sec. VIII H, I calculate (M) and (J) for a slowly- 
moving binary system consisting of bodies of comparable 
masses (one being the black hole); and in Sec. VIII I, I 
take the mass ratio to be small (M/M cxt <C 1) but allow 
the black hole to move rapidly in the strong gravitational 
field of the external body. 

In Sec. IX, I calculate (M) and (J) for a Kerr black hole 
moving in an external universe. I again work to leading 
order in a SH/SM approximation, but the statement of 
the approximation must now be refined to M/TZ <C x, 
where x = a/M = J/M 2 is the dimensionless rotational 
parameter of the black hole. The results, which were 
obtained previously by D'Eath [l^, are 



(M) = 0(M b /H°) 



(1.13) 



and 



(j) = 



3(4+17x 2 )(£ 2 + £ 2 



15 X 2 (E 3 + B 3 



(1.14) 



where E x = £ a p£ afi , E 2 = E^E^ , E 3 = 
(E af} s a s^) 2 , and B 1 = B aP B a ^ B 2 = B a ^B a ^ , 
B3 = (B a ps a s 13 ) 2 . The leading-order calculations car- 
ried out in Sec. IX are not sufficient to determine (M), 
but they indicate that (J) scales as M 5 /TZ 4 . This result 
can also be understood in terms of a coupling between 
the tidal fields and the hole's induced mass and current 
quadrupole moments; here, as I show in Sec. IX E, the 
relationship between M a p and E a p, and the relationship 
between J a p and B a p, do not involve a time derivative 
(as they do in the case of a Schwarzschild black hole). 
Three illustrative applications of Eq. i|1.14fl are worked 
out: In Sec. IX G, I examine a Kerr black hole in circular 
motion in a slowly-moving binary system; in Sec. IX H, 
I consider instead the case of a small hole in relativistic 
circular motion; and in Sec. IX I, the Kerr black hole is 
placed in a static tidal gravitational field. 

The main results of Sec. VIII, Eqs. (|TTT1) and (|l~T2"j) . 
hold to leading order in M/K < 1, and they reveal that 
for a Schwarzschild black hole, (M) = 0(M 6 /TZ 6 ) and 
(J) = 0(M 6 /TZ 5 ). On the other hand, the main results 
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of Sec. IX, Eqs. (|1.13|) and Q1.14JI . hold to leading order 
in M/1Z <C x, and they reveal that for a Kerr black hole, 
(M) = 0(M 5 /K 5 ) and (j) = 0(M 5 /K 4 ). The scal- 
ings are thus very different, and the condition M/1Z <C \ 
implies that the Schwarzschild results cannot straightfor- 
wardly be obtained from the Kerr results in a limit \ 0. 
These scalings indicate that when a rotating and a non- 
rotating black hole are placed in identical environments, 
the rotating hole will absorb larger quantities of energy 
and angular momentum. The agent responsible for this 
enhanced absorption is evidently the hole's rotation, and 
some insight into this matter is offered in Sec. IX F. 

This concludes the summary of the work presented in 
this article. 



F. Organization of the paper 

The rest of the paper contains the derivations of the 
results summarized previously, as well as additional re- 
sults and details. I begin in Sec. II with a description 
of the event horizon of an unperturbed Kerr black hole, 
and I derive a number of results that will be used in 
later sections of the paper. A description of a general 
evolving horizon is presented in Sec. III. This discussion 
is specialized, in Sec. IV, to the case of a perturbation 
of the Kerr horizon, and I derive the flux formulae of 
Eqs. (Q~23)-(£01. 

These formulae are translated into a practical curva- 
ture formalism in Sec. V; this was summarized in Sec. I 
C. In Sec. VI they are translated into a metric formal- 
ism that applies to both rotating and nonrotating black 
holes; in that section I introduce "preferred" and "radia- 
tion" gauges for the metric perturbations, and I establish 
a connection between the fluxes and Isaacson's effective 
stress-energy tensor for gravitational waves |3lL l32 | . The 
metric formalism is specialized to a Schwarzschild black 
hole in Sec. VII; this was summarized in Sec. I D. 

The small- hole/slow- motion approximation is worked 
out in the last two sections of the paper, first for a 
Schwarzschild black hole immersed in an external uni- 
verse (Sec. VIII), and next for a Kerr black hole (Sec. IX). 
These results were summarized in Sec. I E. 

The Appendix contains a brief summary of the the- 
ory of metric perturbations of a Schwarzschild spacetime. 
The material contained in this Appendix is well known, 
but it is convenient to record the main results there be- 
cause they are referred to on many occasions in the arti- 
cle's main body. 

Throughout this work I use geometrized units in which 
c = G = 1, and I adopt the conventions of Misner, 
Thorne, and Wheeler yjl. 



II. KINEMATICS OF THE KERR HORIZON 

To prepare the way for the discussion of dynamical 
event horizons in the next two sections, in this section 



I cover the kinematics of a stationary event horizon de- 
scribed by the Kerr metric. I shall introduce a para- 
metric description of the horizon's null generators, and 
derive from this an intrinsic description of the horizon. 
Part of this discussion will be devoted to the construc- 
tion of null tetrads on the horizon, and a description of 
the (well-known) algebraic structure of the Weyl tensor. 



A. Kerr metric 

Throughout this work the Kerr metric will be written 
in terms of ingoing Kerr coordinates (v,r,9,tp), so that 
its form will be regular on the event horizon. It is given 
by (see, for example, Box 33.2 of Ref. Q, or Sec. 5.3 of 
Ref. [13) 



ds 2 



a 2 sin 2 I 



2a sin 2 9 drdip + 



dv 2 + 2 dvdr - 



AMar sin 2 9 



£ sin 2 9 



■dip 2 



■? 2 d9 2 , 



dvdip 
(2.1) 



where M is the black-hole mass, J = Ma its angular mo- 
mentum, p 2 = r 2 + a 2 cos 2 9, A = r 2 — 2Mr + a 2 , and S = 
(r 2 +a 2 ) 2 — a 2 A sin 2 9. The transformation from the more 
usual Boyer-Lindquist coordinates (£bl, ^bl, #bl> </>bl) is 
given by v = t Bh + J (r 2 + a 2 )A _1 dr, r = r Bh , 9 = 9 BL , 
and ip = 4>bl + o, J A -1 dr; the Kerr coordinates are 
sometimes denoted (V,r,0,<f>), as is done in Ref. |34j . 
The event horizon is situated at the largest root of A, at 
r = r+ = M + VM 2 - a 2 . 

The Kerr spacetime admits the Killing vectors t a = 
dx a /dv and cp a = dx a /dip. The vector 



k a =t a + fli 



with 



2Mn 



(2.2) 



(2.3) 



is also a Killing vector, and it is null on the event hori- 
zon; it is tangent to the horizon's null generators. The 
quantity iljj is the angular velocity of the black hole. 



B. Parametric description of the horizon 

We wish to introduce a system of coordinates (v,9 A ) on 
the horizon, adopting v as a parameter on the generators, 
and 9 (A = 2,3) as generator labels that stay constant 
as the generators move. Because k a — (1, 0, 0, f2n) in the 
spacetime coordinates (v,r,9,ip), we have that 9 is con- 
stant on each generator, and it can therefore be chosen 
as one of the comoving coordinates. On the other hand, 
dip/dv = fin and tp increases linearly as the generators 
wrap around the event horizon; a suitable choice of co- 
moving coordinate is therefore <p = tp — £Ih_v, which stays 
constant. Our horizon coordinates are therefore 



( v ,e A ) = (v,e,<p = i>-n a v). 



(2.4) 



7 



It is important not to confuse the horizon coordinate <j> 
with the Boyer-Lindquist coordinate </>bl; these are not 
equal. 

The horizon generators can now be described by para- 
metric equations of the form x a — z a (v,9 A ), in which 
z a gives the spacetime-coordinate positions of the gener- 
ators in terms of the intrinsic horizon coordinates. Ex- 
plicitly, the parametric description is v = v, r = r+, 
9 = 9, and ip = <j> + fljiv. The vectors 



dz a 
06 1 



are tangent to the horizon; k a is tangent to each gener- 
ator while e A points in the directions transverse to the 
generators. In the spacetime coordinates (v,r,9,ip) we 
have k a = (l,0,0,O H ) as before, e£ = (0,0,1,0), and 
= (0,0,0,1). Because the coordinates 6 A are comov- 
ing, the transverse vectors e A are Lie transported along 
the generators, and they therefore satisfy 



p ot k/3 _ i,a J3 
e A;0 K — K ;0 e A- 



(2.6) 



They are also Lie transported along one another, so that 

P a „0 _ p a J3 

The basis vectors also satisfy 



k^k — — k^Gj^. 
The only nonvanishing inner products are 

1AB = g a peA e Bi 



(2.7) 



(2. 



and these form the components of the induced metric 
on the horizon. To see this, deduce from Eq. (|2.5(1 
that a displacement on the horizon is described by 
dx a = k a dv + e a A d6 A and calculate ds 2 = g a pdx a dx^ 
for this displacement; use of Eqs. (|2.7(l and (|2.8(l re- 
turns ds 2 = ■jab d6 A dd B , with the interpretation that 
jab is indeed the induced metric. Notice that the hori- 
zon metric is degenerate, and explicitly two-dimensional 
in the comoving coordinates. The nonvanishing compo- 
nents of the horizon metric are jgg = ri_ + a 2 cos 2 6, 



(r 2 + + a 2 ) 2 sin 2 9/(r\ + a 2 cos 2 6) 
^/j = (r+ + a 2 ) sin 6* 



and 



(2.9) 



is the square root of the metric determinant. 

The vector basis on the horizon can be completed with 
another null vector N a that satisfies 



where j AB is the inverse of jab- In the sequel we will 
use the horizon metric and its inverse to lower and raise 
upper-case Latin indices. We will also introduce a two- 
dimensional connection T BC compatible with jab, and 
denote covariant differentiation in this connection with a 
vertical stroke; for example, Jab\c = 0- 



by 



Horizon connections 



The tangential derivatives of the basis vectors are given 



ft . I^'v — flirt, , 



k a . 
el 



(3 C A 



io A k a 



= PAB k a + T c AB e a c 



-B;P^Ai 



(2.13) 
(2.14) 
(2.15) 



where k, uia, Pab, and L^ B are the horizon connections 
The surface gravity n 
is given by 



N a k a .pk° of a Kerr black hole 



M /M 2 " 



2Mr 4 



(2.16) 



and Eq. H2.13fl states that the vector k a satisfies the 
geodesic equation, but that the generator parameter v is 
not an affine parameter. Explicit expressions for loa and 
Pab will not be needed; the two-dimensional connection 
^ab can ea sily be computed from jab ■ 

The 2-vector <f) A = (0, 1) is a Killing vector of the 
horizon's intrinsic geometry, and it therefore satisfies 
Killing's equation, 4>(a\b) = 0. This vector is related 
to the spacetime Killing vector <j) a by the relation <j) a = 



(j) A e% The 2-tensor 



cab 



»A\B 



(2.17) 



will be needed in Sec. VI of the paper. This tensor is anti- 
symmetric by virtue of Killing's equation; its only nonva- 
nishing components are cg$ = (r 2 + a 2 ) 3 sin 6 cos 9/{r 2 _ + 
a 2 cos 2 9) 2 = —c^e- 

The vector k a introduced in Eq. I|2.5|l is defined on the 
horizon only, but Eq. (|2.2I) provides an extension away 
from the horizon. The extended vector field is null on 
the horizon only, but it is everywhere a Killing vector; it 
satisfies 



N a N a = = N a e% N a k a = -l. 



(2.10) 



fca;/3 = -K(k a Nf3 - N a kp) + k a UJf3 - LO a k(3 (2.18) 



These conditions determine the vector uniquely, and we 
find 



»r , n , a 2 sin 2 9 , 
N a dx = -dv + — g- dr. 

2{r+ + a z ) 



(2.11) 



The four basis vectors give us completeness relations for 
the inverse metric evaluated on the event horizon, 



on the horizon, where oj a = LO A e° 



D. Null tetrads 



The transverse vectors can be combined into com- 
plex vectors e Q = e A e A that satisfy 



(2.12) 



6q6 







6qC 



6a 6 



(2.19) 
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with an overbar indicating complex conjugation. In 
terms of the complex coefficients e A , these relations read 



1AB& e 



A - B -0 = ~f AB e A e B , lAB e A e B = l, (2.20) 



and these produce the completeness relations r y AB = 
e A e B + e A e B . Substituting this into Eq. I|2.12|l yields 

g a P = -k a N p -N a k p + e a e p + e a e . (2.21) 

A particular (and traditional) choice of coefficients e A 
that achieves these properties is 



1 



lacost 



la cos ( 



V2(rl 



a 2 ) sin# 
(2.22) 

Note that the inversion formula is e A = e A e a + e A e a , 
where the upper-case Latin index was lowered with the 
horizon metric "/ab ■ 

The basis (k a , N a , e a , e a ) is a null tetrad on the hori- 
zon, and it can be used to decompose various tensors, as 
is customary in the Newman-Penrose formalism (see, for 
example, the presentation of Ref. [3q|l. This tetrad, how- 
ever, is not adapted to the algebraic structure of C ai ps, 
the Weyl tensor of the Kerr spacetime. With this tetrad 
we would find that the Weyl scalars tpa & n d tpi vanish, 
but that ip3 and ip4 do not. (These quantities will be 
introduced below.) This is remedied by a null rotation 
(a rotation of class I in the language of Chandrasekhar 
|35j) to a new tetrad (k a ,n a ,m a ,rh a ) given by 



n a = N a + \A\ 2 k a + Ae a + Ae a , m 
where 

. i , „ r + — iacost 
A = -=osmt 
y/2 



Ak a , 
(2.23) 



(2.24) 



In the spacetime coordinates (v,r,9,ip), the components 
of the new transverse vectors are 



1 



la cos t 



rn 



y/2', 



ia sin#, 0, 1, ■ 



smf 



(2.25) 



Notice that this rotation leaves the vector k a unchanged. 

The tetrad (k a ,n a ,m a , m Q ) introduced here is the 
Hartle- Hawking null tetrad of the Kerr spacetime |26j (as 
it was defined by Teukolsky and Teukolsky and Press 
PJ) restricted to the event horizon. The Hartle-Hawking 
tetrad is related to the more standard Kinnersley null 
tetrad [3(| by a rescaling of the vectors k a and n a (a 
rotation of class III in the language of Chandrasekhar 

My- 



fc Q (HH) = 



A 



2(r 2 



-k a (K) 



and 



n Q (HH) 



2(r 2 +a 2 ) 



n"(K). 



(2.26) 



(2.27) 



Notice that while the Hartle-Hawking tetrad is de- 
fined globally in the Kerr spacetime, the tetrad 
(k a , n a , m a , rh a ) is defined on the horizon only; the ex- 
tension of k a to k a (HH) is different from the extension of 
Eq. (|2.2I) . Notice also that whereas the Hartle-Hawking 
tetrad is well behaved on the event horizon, the Kinner- 
sley tetrad is not. 

In the tetrad (k a ,n a ,m a , fh a ) the horizon Weyl scalars 
are defined by (see, for example, Ref. |35jp 



i>o = 


-C aif5S k a m~<kPm s , 


(2.28) 


rh = 


-C^ps^n^k^m 5 , 


(2.29) 


^2 = 


~C ai 0&k a m 1 rhPn h \ 


(2.30) 


V>3 = 


—C ai f3sk a n 1 fh ,3 n 5 , 


(2.31) 


lp4 = 




(2.32) 



where C a ^ps is the Weyl tensor of the Kerr spacetime 
evaluated on the event horizon. We have that ipo = ipi = 
ip3 = ipi — 0, a property that is true globally with the 
Hartle-Hawking or Kinnersley tetrads, and 



f/'2 



M r + (r 2 + - 3a 2 cos 2 9) iMa cos 6»(3r 2 



in + ■ 



(r 2 + a 2 cos 2 9) s 



(2.33) 



is the only nonvanishing horizon Weyl scalar. 

E. Weyl identities 

For future reference I record here a number of identities 
satisfied by the Weyl tensor of the Kerr spacetime: 



Caypgk^ e^k^ 
C ai ps(e A eg + e%e A )k 7 



0, (2.34) 
-2Re(^ 2 ) 1A Bks, (2.35) 
C^psk^k 5 = -2Rc(^ 2 )k a kp. (2.36) 

As a consequence of Eqs. (|2.34|) . I|2.35|l . and 12.2311 we 

also obtain 

C^ps^m^k 5 = = C ai p S m a m l3 k^. (2.37) 
These identities are formulated on the horizon only. 

III. DYNAMICS OF AN EVOLVING HORIZON 

In this section I generalize the preceding discussion to 
a nonstationary event horizon. The presentation is pat- 
terned after Price and Thorne [2?! an d Chapter VI of the 
book by Thorne, Price, and Macdonald |28j . 

A. Comoving coordinates and vector basis 

As in Sec. II we take (v, 9 A ) as our system of intrinsic 
coordinates on the horizon, with v now promoted to an 
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arbitrary parameter on the null generators, and 9 A still 
denoting constant generator labels. The horizon is still 
described by the parametric equations x a = z a (v,9 A ), 
but the coordinate positions of the dynamical horizon 
may be displaced with respect to those of a stationary 
Kerr horizon. 

The vectors k a = dz a /dv and e\ — dz a /d9 A form a 
partial basis on the horizon; as before k a is tangent to 
the generators, e a A is transverse to them, and k a k a = 
k a e\ — 0. The nonvanishing inner products 



1ab(v,9 a ) = g, 



a 
a/3 e A e B 



(3.1) 



give the components of the induced metric on the horizon. 
The basis is completed by another null vector N a , which 
is orthogonal to e A and normalized by the condition 
N a k a = — 1. The completeness relations of Eq. (12.12(1 
still apply. 

The vectors k a and e A are all Lie transported along one 



another, so that e A . k^ = k a ^e A and e A .0e B = e%.^e. A . 



B. Generator kinematics 

The tangent vector k a satisfies the geodesic equation 
in its generalized form 



(3.2) 



where n is the evolving surface gravity of the event hori- 
zon, defined with respect to our choice of parameteriza- 
tion v. 

The transverse derivatives of the tangent vector can be 
decomposed as 



;P e A 



u A k a + B 



Bn 



-A;/3 



k\ 



(3.3) 



for some 2-vector u>a(v, a ) and 2-tensor Bab(v, 9 a ); this 
generalizes Eq. (|2.14|l . It is easy to show that the right- 
hand side of Eq. 1)3.3(1 cannot include a term proportional 
to N a ; this follows from the fact that k a is a null vector 
field. And it can be established that Bab is a symmetric 
tensor, because the congruence of null geodesies to which 
k a is tangent is necessarily hypersurface orthogonal (see, 
for example, Sec. 9.2 of Ref. @, or Sec. 2.4 of Ref. 
The tensor Bab = k a -0e A e B can be decomposed into its 
irreducible parts, 



Bab = -jQlAB 



0~AB, 



(3.4) 



thereby defining the expansion scalar = 7 Bab and 
the shear tensor gab — Bab — \®1ab- Notice that the 
expansion is the trace of Bab, while the shear is the 
tracefree part of this tensor. 



C. Generator dynamics 

Evolution equations can easily be derived for jab, ©j 
and o AB, and these will form the basis of the discussion 
of perturbed event horizons in the next section. 



Starting with the identity d"fAB/dv = ((; Q ^e^e^) ;7 fc 7 
and using Eqs. I|3.3I) and (|3.4|l . we quickly arrive at an 
evolution equation for the horizon metric, 



97. 



AB 



dv 



QjAB + ZCTAB- 



(3.5) 



From this it follows that dj AB /dv = -Oj AB -2a AB and 

1 5V7 







V7 dv 



(3.6) 



where 7 is the metric determinant. 

To derive evolution equations for the expansion and 
shear we follow the usual route that leads to Raychaud- 
huri's equation (see, for example, Sec. 9.2 of Ref. |29| . 
or Sec. 2.4 of Ref. |3(]]). Starting with the identity 
9Bab I dv = (/c Q;(3 e^e|) ;7 /cT and using Eqs. (E3, lpO|l . 
as well as Ricci's identity, we arrive at 8Bab/9v — 
kBab + B A C Bcb — R a ~ f i3se'\k' 1 e%k & . Taking the trace of 
this equation, using the fact that a symmetric-tracefree 
tensor automatically satisfies o- A G a c B = ^{o-cdO' CD )S a b , 
produces Raychaudhuri's equation, 



dv 2 



87175, 



(3.7) 



where p = (R a/} /8n)k a k' :! = T^k -^ after using the 
Einstein field equations. The tracefree part of the equa- 
tion reduces to 



d° A B 
dv 



= ( K -G)a A B -C A , 



where 



uS 



Cab = C arj pse°Ak~ 'e B k 



(3.8) 



(3.9) 



are tangential components of the Weyl tensor. In these 
equations all upper-case Latin indices are lowered and 
raised with jab and "f AB , respectively; the horizon met- 
ric evolves according to Eq. I|3.5(l . 

The area of any cross section v — constant of the event 
horizon is given by A(v) — § ^jd 6. Assuming that 
the number of generators stays constant as the horizon 
evolves (that is, assuming that no new generator joins 
the horizon at a caustic), a change of area occurs when 
7, the metric determinant, varies with time. Equation 
(E3 yields 



^= <peds, 

dv 



(3.10) 



where dS = ^/jd 2 9 is an element of surface area on the 
horizon cross sections. 

The equations derived in this section are all exact, and 
they apply to an event horizon that evolves dynamically. 
[The assumption made in the derivation of Eq. I|3.10|l . 
that the number of generators must be conserved during 
the horizon's evolution, represents a serious restriction. 
In the perturbative context to be described in the next 
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paragraph, however, this limitation is lifted because the 
formation of a caustic is necessarily associated with a 
large, nonperturbative value of O.] The choice of param- 
eter v and generator labels 9 A is completely arbitrary, 
and the quantities K, jab, 0, and gab all refer to this 
choice. For a stationary Kerr black hole, v is chosen so 
that k a is given by Eq. i|2.2[l . and the generator labels 
of Eq. I|2.4[) are adopted. This means that k is given by 
Eq. (|2.16(l . jab by the expressions listed near Eq. 12.9|l . 
and that O = = (Jab, as can be seen by comparing 
Eq. (QUI with Eq. Equations (E3 and (231 are 

then consistent if p = 0, which follows because the Kerr 
metric is a vacuum solution to the Einstein field equa- 
tions, and Cab = 0, which follows from Eq. l|2.35|l . 

In the next section we will apply these equations to 
situations in which the horizon is very close to being sta- 
tionary, so that it can be described as a slightly perturbed 
version of the Kerr horizon. The horizon coordinates 
(v,9 A ) will then be chosen to be "close to" the Kerr co- 
ordinates, and we will see that the ambiguities associated 
with this choice never explicitly enter the discussion; our 
final expressions will be gauge invariant. This implies 
that the vectors k a , N a , e a A on the perturbed horizon 
will be perturbed versions of the Kerr basis, and that all 
derived quantities, such as k, jab, 0, and gab, will be 
perturbed versions of the corresponding Kerr quantities. 



IV. PERTURBED HORIZON 

I now specialize the formalism of the preceding section 
to event horizons that are slightly nonstationary, those 
that can be considered to be perturbed versions of the 
Kerr horizon. I shall assume that the perturbation is 
caused entirely by gravitational radiation, and no matter 
will be allowed to cross the event horizon. The pertur- 
bation formalism described here is adapted from Price 
and Thorne [53] and Chapter VI of the book by Thorne, 
Price, and Macdonald |2Sj: these methods go back to the 
pioneering work of Hawking and Hartle 26] . 



A. Perturbation equations 

Changing our notation with respect to the previous 
section, the perturbed values for the horizon metric, sur- 
face gravity, expansion, shear, and Weyl tensor will now 
be denoted jab, k, O, gab, and Cab, respectively; these 
quantities were all introduced in Sec. III. The unper- 
turbed (Kerr) values will be denoted without a decorat- 
ing caret; for example = is the background expan- 
sion scalar, gab — the background shear tensor, and 
Cab — the background Weyl tensor. The only non- 
vanishing background quantities are the metric jab and 
surface gravity k; these were introduced in Sec. II. 

The horizon perturbation is driven by the Weyl tensor 
Cab, which we imagine to be a quantity of the first order 



in an expansion parameter A; we write this as 
C A B = \C A B + 0{\ 2 ). 



(4.1) 



At the end of the calculation we will set A = 1 by absorb- 
ing it into the definition of the perturbations. Equation 
(|3.8|) indicates that the Weyl tensor drives a first-order 
perturbation in the shear, and we have 



G A B =\G 1 A B +0{\ 2 ). 



(4.2) 



Equation (|3.7() . on the other hand, shows that it is the 
square of the shear tensor that is driving a perturbation 
in the expansion (recall that we have set p = 0), and we 
must therefore have 



= A 2 2 + 0(A 3 



(4.3) 



(A more careful treatment that incorporates a first-order 
term would eventually lead to Oi = and therefore back 
to this assertion.) Finally, Eq. (|3.5|l shows that the shear 
produces a first-order perturbation in the metric, 



JAB = JAB + ^JAB +0{\ 2 ), 



(4.4) 



and these considerations lead to the statement that the 
perturbed surface gravity will differ from its Kerr value 
by a first-order quantity: k = n + O(A). 

Substituting the expansions of Eqs. I|4.1|) - (|4.4|) into 
Eqs. I|3.5|l . I|3.7|l . and (|3.8() gives us the following set of 
perturbation equations: 



Ha 



AB 



dv 
502 
dv 

dfjl AB 

dv 



= 2cr AB, 



K&2 — a\ 



AB 
AB U \ ) 



— kg\ b — C AB . 



(4.5) 
(4.6) 
(4.7) 



In these equations, k is the surface gravity of the Kerr 
black hole, and the generator parameter v can be identi- 
fied with the advanced-time coordinate of the Kerr space- 
time (as introduced in Sec. II); the first-order deviations 
of these quantities with respect to the Kerr values do not 
enter the perturbation equations. It should be noted that 
upper-case Latin indices can now be manipulated with 
jab , the background horizon metric, and its inverse. 

According to Eq. (|3.10|) . a growth in the horizon area 
is driven by the expansion scalar, and we therefore have 
A = A 2 i 2 + 0(A 3 ), with 



A 2 = ^ = ./ 0,,/.S'. 



(4.8) 



where dS = ^Jjd 2 9 is an element of surface area on the 
cross sections of the unperturbed (Kerr) horizon. 

B. Integration of the perturbation equations 

We imagine a horizon that starts in an initial Kerr 
state, is perturbed for some time by an external process, 
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and then eventually returns to another Kerr state. The 
expansion and shear vanish in the initial state, and they 
must return to zero after the external process has ended; 
this requires the imposition of teleological boundary con- 
ditions (as opposed to retarded boundary conditions; see, 
for example, Sec. VI C 6 of Ref. 28] ) on the solutions to 
the perturbation equations. 

The teleological solution to the equation (d/dv — k)i/j — 
-f(v) is if>(v) = e' K{v '- v) f(v')dv'. It shows that 
ip(v) depends on the future behavior of the driving force, 
but that ip{ v ) S oes to zero after the driving force is 
switched off; the causal solution would depend only on 
the past behavior of the driving force, but it would grow 
exponentially after the force is switched off. If the driving 
force varies very slowly over a time comparable to 1/k, 
then the teleological solution reduces to the local expres- 
sion ip(v) = k _1 /(w), to a fractional accuracy of order 
(kt) _1 , where t ~ /// is the time scale over which the 
driving force varies. The local expression can be simply 
obtained by noting that in this limit, d/dv <C K and the 
differential term can be neglected in the differential equa- 
tion; a more careful derivation starts with the teleological 
solution and employs integration by parts. 

It is not permissible to neglect the differential term in 
Eq. H4.7f) . and one must write down a proper teleological 
solution to this equation. To see this, suppose that the 
black hole is a member of a binary system, and that it 
moves in the field of an external body with an angular 
velocity f2. As seen in the rotating frame of the gen- 



erators, the Weyl tensor behaves as C 



where 



uj = f2 — Qu is the relative angular velocity between the 
external field and the generators. Thus, unless the exter- 
nal field is nearly corotating with the black hole, u) will 
be of order J1h, which is itself of order k, and the Weyl 
tensor will not vary slowly. 

The exact solution to Eq. 14. 71) is 



^ab(v,0 A ) 



e -W-) C \ B { v ',6 A )dv', (4.9) 



and this can be substituted into the right-hand side of 
Eq. I|4.6|) . Here we shall allow ourselves some simplifi- 
cation. For a Weyl tensor that behaves schematically as 
C ~ e~ luJV , the shear tensor will go as a ~ e~ %u)v / (n+iuj) , 
and it will vary as rapidly as C . The square of the shear 
tensor, however, will contain a piece that oscillates at 
twice the frequency u>, and a piece that stays constant. 
The force driving the expansion therefore contains both 
a slowly- varying piece and a rapidly- varying piece. In a 
generic situation we expect that o~ AB a AB can always be 
decomposed into such slowly- varying and rapidly- varying 
pieces, and we isolate the slowly- varying component by 
averaging over a time scale that is long compared with 

K ~ 1[ (°~AB a i B )si°™ = ( (J AB cr i B )- If we now agree to 
follow only the slow evolution of the expansion scalar, 
and to ignore its rapid fluctuations around the slowly- 
evolving mean, then 



1 



<e 2 ) = >W B > 



(4.10) 



is an adequate approximate solution to Eq. I|4.6|) . The 
meaning of the averaging sign should be clear: The ex- 
pansion scalar, and the square of the shear tensor, are 
averaged over a time t that is long compared with ft -1 , 
the black-hole time scale. The time scale r is identified 
with a characteristic time associated with the growth of 
the black-hole area, r ~ (A) / (A) . We note that it is a 
requirement of the perturbative treatment that kt ^> 1; 
the simplification of Eq. i|4.10[) therefore represents no 
significant loss of generality, other than a coarse-graining 
over short time scales. 

Equation l|4.9|l can also be inserted into Eq. 14.511 to 
calculate the metric perturbation, which is given by 



' AB 



(v',9 A )dv'. 



(4.11) 



After altering the order of the integrations and perform- 
ing one of the integrals, we obtain 



1ab(v,0 A ) = 



c\ B {v',e A )dv' 

e-< v '-^C\ B {v',e A )dv' 



■ (4.12) 



this result is exact, and it does not involve a coarse- 
graining over short time scales. 

The averaged rate of change of the horizon area can 
be calculated on the basis of Eqs. (|4.8|) and (|4.10|l . The 
result is 



(A 2 



I fi^AB*?*) dS, (4.13) 



where an overdot indicates differentiation with respect 
to v. This result can be expressed in terms of the Weyl 
tensor by means of Eq. Q4.9f) . 

In the remainder of the paper we will denote C\ B sim- 
ply as Cab, o~\b as °~ab, ©2 as G, and A 2 as A; because 
these quantities all vanish for a stationary horizon, this 
change of notation will not generate ambiguities. But to 
avoid ambiguities we will continue to denote the metric 
perturbation as j\ B . 



C. Fluxes of mass and angular momentum 

We shall now derive expressions for the averaged rates 
of change of the black-hole mass M and angular momen- 
tum J, using Eq. (|4.13|l as our main input. 

In the case of a horizon perturbed by a matter field it 
can be shown (see, for example, Sec. 6.4.2 of Ref. |37|1 
that these rates are related by (k/8it)A = M — QhJ, 
which is a statement of the first law of black-hole mechan- 
ics (see, for example, Chapter 12 of Ref. [2j|, or Chapter 
5 of Ref. [13 )• In the present case of a horizon perturbed 
by a purely gravitational perturbation, we shall assume 
that this relation holds on the average, so that 



8tt 



(A) = (m)-Qh(J). 



(4.14) 
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The averaging introduced here is the same that was in- 
volved in Eq. (|4.1U|I . If we imagine that the horizon 
evolves from an initial Kerr state to a final Kerr state in 
a time Av <~ r, then a precise statement of Eq. (|4.14|) is 
(k/8tt)(SA)/(Av) = (SM)/(Av) - tt K (5J)/(Av), where 
5 A, 6M, and 5 J are the total accumulated changes in 
the black-hole parameters. Because these changes relate 
two stationary black-hole states, we have here the usual 
statement of the first law divided by the time interval 
Aw. 

In the case of a horizon perturbed by a matter field 
it can also be shown (see, for example, Sec. 6.4.2 of 
Ref. [33) that if the matter field is decomposed into 
modes proportional to e~ wv e im ^ , where a; is a contin- 
uous frequency and to an integer, then each mode con- 
tribution to the averaged rates is such that 

(M)m,o, = —(j) m , u . (4.15) 

TO 

This mode decomposition is motivated by the symmetries 
of the background Kerr spacetime, and Eq. Ij4.15|l states 
that a mode labeled by (to, lu) carries across the horizon 
a quantity of energy proportional to u> and a quantity of 
angular momentum proportional to to. This statement is 
easily understood on the basis of a quantum picture, but 
it holds for classical matter fields as well. We shall as- 
sume that Eq. (|4.15l) is not restricted to matter fields, but 
that it holds also for gravitational perturbations. Such 
an assumption was made previously by Teukolsky and 
Press |l| in their pioneering study of horizon fluxes. 
Equations (|4.14() and (|4.15|) imply 

(M)m,w = (4.16) 
k 8n 

(j)m^ = ^^-{A) m , ui (4.17) 
k 8n 

where 

k = Lu - to^h (4.18) 

and {A) mtUJ is the mode contribution to the averaged rate 
of change of the horizon area. These equations will allow 
us to turn Eq. I|4.13|) into useful expressions for (M) and 
(j). 

In the spacetime coordinates (v,r,9,ip), the mode de- 
composition of the metric perturbation is given by 

^ab = E / ^AB(r, 0)e~^e^. (4.19) 

In the horizon coordinates (i>, 9, 4>), where <p = ip — rt^v, 
we have instead 

Iab = E/ ^2B(r + ,6)e-^e™t (4.20) 

m 

where we have set r = ?*+, choosing the coordinate po- 
sition of the perturbed horizon to coincide with the po- 
sition of the Kerr horizon. (That this choice can always 



be made is proved in Sec. VI A.) Substituting Eq. (|4.20() 
into Eq. 14. 5|) we obtain 

= ^ E / dui-ik^e-^e^t (4.21) 

m 

the mode decomposition of the shear tensor. 

We now insert Eq. (|4.21[) into Eq. 1|4.13[) . but we do 
not decompose a AB into modes. This gives 

£<4 = £/ d - iS fi^i&e-^e**) dS, 

m 

and from this we can read off each mode contribution to 
(k/8%){A). According to Eqs. lj4"Toj) and (|4"T7|) . then, 
we have 

(M) = E / Ua^^e-^e^f) dS 

and 

( J ) = E / ^TeV Uo AB ^e-^e^)ds. 

m 

The metric perturbation can now be reconstructed from 
its mode decomposition. It is easy to show that the fac- 
tor of — ioj is generated by applying the differential oper- 
ator d/dv — Qnd/d<f> to 7^#(v, 9, 4>)', this can be written 
as a Lie derivative in the direction of t a = k a — ttn<j) a , 
the timelike Killing vector of the background Kerr space- 
time. Similarly, the factor of im is generated by acting 
with d/d(fi, which is a Lie derivative in the direction of 
cf) a , the rotational Killing vector of the background Kerr 
spacetime. 

The final expressions are 

{M) = ^j(a AB £ tl \ B )dS, (4.22) 

(j) = -^<f(v AB £<M 1 AB)dS, (4.23) 

and 

^{A) = ^j>(o AB £ kl \ B )dS, (4.24) 

where we have used aAB — \£kl\s- As we nave seen, 
the Lie derivatives acting on the metric perturbations 
refer to specific directions in the background Kerr space- 
time. In the horizon coordinates [v, 9, </>), the Lie deriva- 
tives take the explicit form 

MIL <4 ' 25) 

and £t — £k ^ On the other hand, in the space- 

time coordinates (v,r,9,ifj) they take the form 

£t=(-Pl , £4>=(St) . (4-26) 
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and £ k — £ t + SIh^V We also note that in Eqs. I|4.22f] - 
(|4.24f) . the surface element dS = ^d 2 d refers to the 
metric 7,4s of the unperturbed Kerr horizon. 

Equations H4.22fl ~ H4.24fl are an excellent starting point 
for the development of a practical formalism to calculate 
the horizon fluxes, a topic we shall turn to in the next 
three sections. These equations are not new: they appear 
in Sec. VI C 11 of the book by Thorne, Price, and Mac- 
donald [2^. The derivation presented here, however, is 
substantially different from theirs, and it incorporates the 
Teukolsky-Press assumption of Eq. if4.15fl : this assump- 
tion was not part of the original derivation, and their 
route from Eq. (|4.13fi to Eqs. (|4.22l) . 14.2311 is not as di- 
rect. It should be clear that while the present derivation 
relies on a mode decomposition of the metric perturba- 
tion, the final expressions involve differential (as opposed 
to algebraic) operations and are independent of the de- 
composition. 

D. Rigid rotation 

The perturbed black hole is part of a system in rigid 
rotation when the vector 

C = t a + VLcj) a (4.27) 

is a Killing vector of both the background Kerr spacetime 
and the perturbed spacetime; here Q is a constant angu- 
lar velocity. An example of a system in rigid rotation 
is when the black hole is a member of a binary system, 
moving with a uniform angular velocity in the field of the 
external body. The fact that is a Killing vector of the 
perturbed spacetime means that 

£ a \ B = 0. (4.28) 

If the metric perturbation is expressed in the spacetime 
coordinates (v,r,9,tp), then Eq. (|4.28fl implies that its 
dependence on v and tp is through the combination ip—Qv 
only; in a reference frame that is rotating at an angular 
velocity £1 with respect to the original inertial frame, the 
perturbation appears to be stationary, and the system is 
indeed in rigid rotation. 

Equations (23 and l|4"27|) imply k a = £ a + (n H -ty<f> a 
and t a = - fl(j) a , and by virtue of Eq. (|4~2g|) the Lie 
derivatives of j\b in the directions of k a and t a can be 
expressed in terms of a derivative along cf> a , 

£ klAB = (^H - ti)£<t>lAB 

and 

£tlAB = -Q £ <plAB- 

Substituting this into Eqs. H4.22fi - (|4.24fi . and using also 
a AB = |(0 H - Q)£^f B , we obtain 

(M) = 0(Q - H )/C, (4.29) 
( j) = (0 - H )/C, (4.30) 
-^(i) = (fi-fi H ) 2 /C, (4.31) 



where 

K=T^j({£^ B ){£^\ B ))dS: (4.32) 

we recall that jf~ B is obtained from r y 1 AB by raising indices 
with 7 , the inverse background metric. 

Equation If4.30[) indicates that the black hole's angular 
momentum will increase when > f2u, that is, when the 
external rotation is faster than the rotation of the gen- 
erators; otherwise the angular momentum will decrease. 
Equation If4.29fl indicates that (M) = fl (j) when the 
black hole is in rigid rotation; the sign of (M) is tied to 
the sign of (J) and the sign of the angular velocity (which 
is denned relative to fin). Finally, Eq. (f4.31fl shows that 
the black-hole area will always increase, as is dictated by 
Hawking's area theorem (see, for example, Sec. 12.2 of 
Ref. (23). These equations also appear in Sec. VII B 1 
of the book by Thorne, Price, and Macdonald 

V. CURVATURE FORMALISM 

In this section I translate the flux formulae of 
Eqs. 14.22fi 14.2411 into a more practical language that 
involves curvature variables. The most important vari- 
able in this formalism is the Weyl scalar ipo , which can 
be obtained by solving Tcukolsky's differential equation 


A. Relation between Cab and tfjo 

Our first task is to express Cab, the Weyl tensor of 
the perturbed horizon, in terms of the more practical 
curvature variable ipo ■ The calculation is straightforward 
but somewhat lengthy; it requires a number of steps. 

The Weyl tensor of the perturbed spacetime is C ai ps — 
Ca-yps + ^C^ 1/3S + 0(A 2 ), and the basis vectors of the 

perturbed horizon are k = k a + Xkf + 0(X 2 ) and = 
e A + Xe^A + 0(X 2 ). The Weyl tensor of the perturbed 
horizon is then defined by 

Cab = Ca^pseA^egk 6 , (5.1) 

which is the same equation as Eq. I|3.9fl . By virtue of 
Eq. If2.35fl we have that Cab = for the Kerr horizon, 
and Cab — ^C\ B + 0(X 2 ). To comply with preceding 
usage we shall now omit the label "1" on the Weyl tensor 
and set A = 1 after the expansion in powers of A has been 
carried out. 

Direct evaluation of Cab from the preceding informa- 
tion gives 

Cab = C^^eAk? 'e^k 5 + C ai ps(e A k? ' e B k{ 

~\~e Ak~^e B k ~\~ eAk^e-^ B k ~{~ e-^A^^ ^b^ )■ 

We will now simplify this expression, and show that Cab 
can be expressed solely in terms of Cl, 1/3S , the perturba- 
tion of the Weyl tensor. 
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The main source of simplification comes from the al- 
gebraic properties of the unperturbed Weyl tensor: The 
last two terms within the brackets vanish by virtue of 
Eq. I|2.34[l . and the first two can be rewritten with the 
help of Eq. fHSty . This gives 



C 



AD 



Cl^s^e^k 5 - 2Re(^ 2 ) 7AB fc Q fc«. 



Next we decompose e\ in the basis of complex vectors 
(e",e Q ) introduced in Sec. II D; the relations are e\ — 
e.A& a + eAe a with the coefficients ca = ^ab& B obtained 
from Eq. TTTfy . This yields 

Cl lf3S eW4k 5 = e A e B Ci 7(3S e a k^e f:! k s 

+ eAeBCl^k^e^k 8 

+ (e A e B + e A e B )^^e a Fe^ 5 

for the first term in the previous expression for Cab- 
The factor within brackets is 7ab, and by virtue of the 
symmetries of the Weyl tensor and the completeness re- 
lations of Eq. J22TJ), we also have Cl 7/3S e a k^e^k s = 

\g al3 C^pgk 1 k s , where g a<3 is the inverse of the Kerr 
metric. But because the perturbed Weyl tensor must 
be traceless in the perturbed metric, = g a ^C ai ps = 
{g al3 -\h a P){C ail 3 S + \Cl ll3S ), where \h aj3 = g a p-g a p is 
the metric perturbation. We therefore have 5 Q/3 C^ 7 ^ (5 = 
h a ^C a -fPS, and gathering these results we obtain 

Cl l(i& e a A We B k 8 = eAeBC^e^e^k 5 

+ eAeBCi^k^e^k 5 

+ \lABh ap C ail)s k^k 5 . 

We may simplify this further by using Eq. Ij2.36|l , and we 
now have 

Cab = e A e B C 1 a ^ s e a k-'e l3 k s + eAeBCi^^k^^k 6 



2Re^ 2 hAB ( ^h a0 k a kP + k a k^ 



The last two terms vanish by virtue of Eq. I|2.37[l . and 
we have —ipo = C„ jggk a m" f k^m s . We now express m a 
in terms e a and k a , as in Eq. I|2.23[l . This yields 



(5.4) 



as our final expression for the Weyl scalar. 

The relation between Cab and "00 is obtained by sub- 
stituting Eq. (|5.4|) into Eq. (|5.2|) . We write this as 



Cab = e A eB^ + e^es*, 



where 



*(«,^)e-4(HH)e- 



A 2 



4(r 2 + a 2 )' 



(5.5) 



rVo(K) (5.6) 



is our main curvature variable. We choose, for conve- 
nience, to absorb a minus sign into the definition of ^ . 
The first equality in Eq. I|5.6|l states that apart from this 
minus sign, '5 is the Weyl scalar ipo as defined with the 
Hartle-Hawking tetrad [2, , evaluated on the horizon 
and expressed in terms of the horizon coordinates. The 
second equality gives the relationship between ^ and the 
Weyl scalar as defined with the Kinnersley tetrad 0, H|| , 
evaluated in the limit r —> r + in which tJjq (K) diverges as 
A -2 . 

The Weyl scalars in either choice of tetrad, and there- 
fore 'J', can be obtained by solving the Teukolsky equa- 
tion. Because VP is a scalar that vanishes in the Kerr 
spacetime, this quantity is gauge invariant. These two 
properties make ^>(v,9 A ) an especially useful choice of 
variable to describe the horizon perturbation. 



B. Fluxes 

We obtain the shear tensor by inserting the Weyl ten- 
sor of Eq. (|5.5() into Eq. I|4.9|) . Because the transverse 
vectors depend on 9 A only, we obtain 



In the last step we recognize that the vector k must be 
null in the metric g a /3, so that = (g a p + Xh a p){k a + 
AfcfX^+Afcf ) = X(2k a k^ + h al3 k a k 13 ). We finally arrive wne re 
at the expression 

Cab = e A e B Ci^ s e a k^e p k s + e A e B Cl l!3 ,e a k?^k\ *+(«> J 

which involves only the perturbed Weyl tensor and the 
Kerr basis vectors. 

The Weyl scalar of the perturbed spacetime is defined 
as in Eq. (|2~2^|) . 



<jab{v, a ) = e A e B $+ + e A e B ^-\ 



(5.7) 



-^'- v )^(v',9 A )dv' (5. 



ipo = C ai ^ik a m 1 k i3 w! > 



(5.3) 



is the future integral of the Weyl scalar iff, weighted by 
the exponential factor e~ K ( v so that only the near 
future contributes significantly. 

The metric perturbation is obtained by substituting 
Eq. I|5.5|l into Eq. I|4.12|l . Here we obtain 



Expansion in powers of A gives -0o = Af/>o + 0(A 2 ). After 
dropping the label "1" and setting A = 1, we obtain 

-0 O = C 1 aif}S k a fn'k p m s + 2C aiP sk a m?k^m{ 

+ 2C a7l 3sk a m~<kfm 5 . 



lAB( v >9 A ) = ^(eAeB$ + e A eB®), (5.9) 



where 



$(v,6 A ) = $-(v,6 A ) + $ + (v,0 A ), (5.10) 
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with 



ty{v',6 A )dv' (5.11) 



representing the past integral of the Weyl scalar 
(weighted uniformly). 

Equations (|5.7I) and l|5.9|l can now be substituted into 
Eqs. (|4.22() - (|4.24l) to obtain expressions for the fluxes. 
Because the basis vectors &a do not depend on v nor (f>, we 
have that £k&A = £cj>eA = £t£A = 0, and the derivative 
operators act only on $±. Using the properties eA& A = 
and eAe A = 1 of the basis vectors — recall Eq. Ij2.20jl - 
we obtain, for example 

(M) = j>(<S>+£t® + dS. 

Using now Eq. Ij2.9jl , we arrive at the following expression 
for the mass flux: 



(M) 



8lTK 



($ + £ t $ + $ + £ f $) dn, (5.12) 



where d£l = sin 9 d9d<f> is an element of solid angle on the 
unit two-sphere. We obtain similarly 



(J) = - 



8lTK 



j {$+£$$ + $+£$$) d£l, 



for the flux of angular momentum, and 

'<l*+| J ><i!i 



s<4 



r\ + a 2 



4tt 



(5.13) 



(5.14) 



for the rate of increase of the horizon area. 

In Eqs. H5.12jl - lj5.14jl it is understood that the 
integrated-curvature fields <§>±(v,0 A ) are expressed in 
terms of the horizon coordinates, and that the Lie- 
derivative operators take the form given by Eq. Ij4.25jl . At 
a later stage we will remove this remaining dependence 
on the horizon coordinates, and express all quantities in 
terms of the original spacetime coordinates. 



C. Pure mode; comparison with Teukolsky and 
Press 

Suppose that the Weyl scalar of Eq. (j5.6|) has the form 

r+,9, ijj) = *t».«(e) e -*-« e <w* (5.15) 

when expressed in terms of the spacetime coordinates; 
this solution to the Teukolsky equation is then a pure 
mode of frequency u> and azimuthal number m. In terms 
of the horizon coordinates we have 



(5.16) 



we recall that <f> = ip — Qnv, and k = uj — mflu was first 
introduced in Eq. Ij4.18jl . We wish to calculate the rates 



of change of mass, angular momentum, and area for this 
pure mode, and to compare our results with those first 
obtained by Teukolsky and Press pj. 

Substituting Eq. (j5.16jl into Eqs. Ij5.8j) . Ij5.11jl . and 
(|5~TU|) yields 



and 



$ = 



# m > w (60 - t ■ 

\ J ^— tkv ^iraq> 

k + ik 
^> m ' UJ (9) t ■ ± 

V / —ikv imd) 

—Ik~ e e ' 



-ik{n + ik) 



(5.17) 
(5.18) 

(5.19) 



It should be noted that while the integral defining $ + 
converges properly for a pure mode, the integral defining 
$_ diverges in the infinite past; this difficulty is remedied 
by inserting a converging factor inside the integral, for 
example e v l Vl with v\ ^> k~ x , to reflect the fact that 
the pure mode was turned on in the finite (but remote) 
past. 

These expressions can now be substituted into 
Eqs. (j5.12|) Ij5.14|) . When acting on a pure mode, £t 
produces a factor of — ilu, £^, a factor of im, and £ k a 
factor of —ik. A simple computation gives 



(M) 
(J) 



—(A) 



w{r\ + a 2 ) 1 
k(K 2 + k 2 ) Air 
m(r^ + a 2 ) 1 
k(n 2 + k 2 ) 4tt 



J (9)\ dil, (5.20) 



j\y m ^(6)\ 2 dSl, (5.21) 



r 2 + a 2 1 



t 2 + k 2 4ir 



* m ' w (6»)f <m. 



(5.22) 



These relations are compatible with Eqs. (j4.29l) - (j4.31|) 

if we define the angular velocity of the pure mode to 
be fl = ui/m; this follows from the fact that according 
to Eq. Ij5.15jl . '5 depends on ip and v only through the 
combination %p—{uj/m)v, so that the perturbation rotates 
rigidly with an angular velocity uj/m. 

To compare our expressions with those of Teukolsky 
and Press pj we must first reconcile the different nota- 
tions. In their Eq. (4.40), Teukolsky and Press display 
the near- horizon behavior of V'o(K) — this is the Weyl 
scalar as defined with the Kinnersley tetrad — in terms of 
Boyer-Lindquist coordinates (£bl, t, 9, 4>bl)- They have 



^o(K) 



-ikr* 



~ 2 2Sh 



hole ; 



where r* = j(r 2 + a 2 )A~ 1 dr = v - ieL, 2Si m (9) are the 
Teukolsky angular functions, and Ihoic is a normalization 
factor. It is easy to check that after a transformation to 
the well-behaved Kerr coordinates (v,r,9,7p), and after 
the rescaling of Eq. (j5.6jl , this expression becomes 



hole 



4(r 2 _ + a 2 ) 2 
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where /3(f) = — a(ri + a 2 ) -1 J(r + r+)(r — r-)^ 1 dr is 
a function that is well behaved at r = r + . Here, r_ = 
M — V M 2 — a 2 denotes the position of the inner horizon. 

The preceding equation relates our definition of 
^ m ^{9) in Eq. I|5.15[) with the quantities introduced by 
Teukolsky and Press. Inserting this into Eqs. (|5.2U[) and 
(|5.21|1 returns 

/d 2 M\ \2Si m (0)Y holc \ 2 UJ 



\ dvdQ I Mv(r\ + a 2 ) 3 k(n 2 + k 2 ) 



and 



\ dvdn I 



\lSi m {9)Y\ 1Q \ e \ 



647r(r^ + a 2 ) 3 k(n 2 + k 2 ) ' 



These are precisely the results obtained by Teukolsky and 
Press |l| and displayed in their Eq. (4.44). Our formal- 
ism (which is based partly on their work) is therefore 
consistent with theirs. 



D. Decomposition into azimuthal modes 

Our aim in this subsection is to derive practical flux 
formulae that are formulated in the time domain, in 
terms of fields expressed in the spacetime coordinates 
(v, r, #/(/>). We shall not, therefore, follow Teukolsky and 



Press [l| and decompose into modes proportional to 
e -!u» e im^ SQ ag j. Q wor k m t ne frequency domain. But 

we will still decompose the Weyl scalar into azimuthal 
modes proportional to e lm ^ , and write 



f(v,r + ,6>,V) = ^* m («,6»)e 



imip 



(5.23) 



This decomposition is motivated by the axial symmetry 
of the Kerr spacetime, which implies that each mode la- 
beled by m will evolve independently Such a decompo- 
sition is therefore likely to be involved in most attempts 
to integrate the Teukolsky equation numerically, in the 
time domain. It will also allow us to remove the re- 
maining dependence of our flux formulae on the horizon 
coordinates (v, 9 A ). It should be noted that the flux for- 
mulae of Eqs. H5.12f> - H5.14|) do not require '3/ to be decom- 
posed into modes; they are therefore ready to be used in 
situations where an azimuthal decomposition is not at- 
tempted. But the implementation of these formulae is 
delicate, because <I>± must be evaluated by integrating 
^(w', r+, 9, i/j) along the horizon generators (integrating 
over dv' keeping <f> = ip — 17hV constant) . Our azimuthal 
decomposition will accomplish this automatically. 

Substituting Eq. I|5.23[) into Eqs. (|5.8|) and I5.11|) gives 
us the azimuthal decomposition of the integrated curva- 
tures, which we express in terms of the horizon coordinate 
4> instead of the spacetime angle ip: 



im<f> 



(5.24) 



where 

$™(v,6>) = e h 
and 

$ m (v,6>) : 



e- ( - K -" nn ^ v 'v m (v',9)dv' (5.25) 



e mnHV '^ m (v',6)dv'. 



(5.26) 



Notice the presence of the oscillating factor e imr2Hl ' 
within the integrals; this comes from the transformation 
between <\> and ip and it reflects the fact that the genera- 
tors wrap around the horizon as v' is integrated forward. 
Notice also that Eqs. I|5.25|l - H5.26|) are now independent 
of <j) or tp, so that it is no longer necessary to specify 
which is to remain constant during integration. 

Equation l|5.24|l can now be substituted into 
Eqs. I|5.12|l - H5.14|) . in which $ = + In the horizon 
coordinates (v , 9, <fi) the operator £k is a partial deriva- 
tive with respect to v, £^ produces a factor of im, and 
£t = £k — ^n£ct>- Simple algebra and integration over 
d<p give 



(M) = 



(j) 



4h' 



E 



2k J <|$+| 2 ) sin 8d9 
imfluj (<§™$ m - $™$ m ) sin 9d9 

,2 i „2 



(5.27) 



4k 



J _ $™|>™) S in 9 d9 7 (5.28) 



and 



(A) = -(rl +a 2 )V /(|$™| 2 )sin^. (5.29) 



8tt 



These are the final form of the flux formulae. Notice that 
these no longer involve the angles <j) and -0, and that all 
fields are expressed in terms of v and 9, coordinates that 
are shared by the spacetime and the horizon. 

The steps required to compute (A/), (J), and (A) 
are therefore these (see also Sec. I C): First, solve the 
Teukolsky equation Q for the functions <3> m (v, 9) defined 
by Eq. I|5.23|l . for all relevant values of m; recall from 
Eq. H5.6fl that ^ is (minus) the Weyl scalar ipo in the 
Hartle-Hawking tetrad, evaluated at r — r + . Second, 
compute the integrals of Eqs. (|5.25|l and (|5.26|l to obtain 
$™ (v, 9). Third, and finally, substitute these values into 
the flux formulae of Eqs. (|5.27l) - (|5.29|) . integrate over d9, 
and sum over m. 



VI. METRIC FORMALISM FOR GENERAL 
BLACK HOLES 

In this section I translate the flux formulae of 
Eqs. 14.22H - H4.24fl into a more practical language that 
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involves metric variables. This translation is most useful 
in the context of a Schwarzschild black hole, for which 
the theory of metric perturbations is well developed; I 
shall consider this specific case in the following section. 
In this section I keep the discussion general, so that it 
applies to both rotating and nonrotating black holes. 

A. Preferred gauge 

We expand the metric of the perturbed black hole as 

g a p = g a p + Xh a p, (6.1) 

where g a p is the metric of the unperturbed spacetime — 
the Kerr metric — and \h a p is the perturbation. (As 
we did previously, we keep A for book-keeping but we 
set it equal to unity at the end of the calculation.) We 
wish first to impose a number of gauge conditions on h a /3, 
which will simplify its relationship with the quantity j AB 
introduced in Sec. IV A. 

Our preferred gauge is a "horizon-locking gauge;" it 
has the property that the coordinate positions of the per- 
turbed horizon are the same as those of the unperturbed 
(Kerr) horizon. As we shall see below, it is always possi- 
ble to make this choice of gauge. In the preferred gauge 
the parametric description of the horizon generators is 
given by 

z a ( v ,e A ) = z a ( v ,e A ), (6.2) 

with z a (v,6 A ) giving the parametric description of the 
Kerr generators. This equality implies 

k a = k a , e a A = e% (6.3) 

so that the perturbation of the tangent vectors is iden- 
tically zero: kf = = ef A in the notation of Sec. V 
A. 

The vector k a must be null, and orthogonal to e A , in 
the metrics g a p and g a p. This observation gives rise to 
the three gauge conditions 

h a pk a k® = = h a pk a e A (preferred gauge). (6.4) 

These equations hold on the horizon only, and we shall 
not need to extend them beyond the horizon. The pre- 
ferred gauge is only partially determined, and the space 
of transformations within the preferred-gauge class is 
large. 

For some purposes it will be convenient to supplement 
the gauge conditions of Eq. i|6.4[l with a fourth condition, 
ha^Ni 3 = 0, so that we have the four conditions 

hapk = (radiation gauge). (6-5) 

These conditions also hold on the horizon only, and again 
we have a large space of gauge transformations within the 
radiation-gauge class. The radiation gauge of Eq. (|6.5|) 



is similar to the one introduced by Chrzanowski (his in- 
going radiation gauge 38] ); but it is distinct because 
Chrzanowski imposes Eq. (|6.5|) as well as g a ^h a p = 
globally in the Kerr spacetime [extending k a away from 
the horizon as fc Q (HH), the first member of the Hartle- 
Hawking tetrad]. The Chrzanowski radiation gauge is 
therefore much more rigidly defined than (and a special 
case of) the radiation gauge of Eq. (|6.5|) . 

Equation l|6.3|l states that three of the basis vectors on 
the horizon are not changed by the perturbation. The 
fourth basis vector, N a , must be orthogonal to e A in the 
perturbed metric, and it must also satisfy g a /3N a k^ — 
— 1 . It is easy to see that these requirements imply N a — 
N a - |A(^ 7 iV^iVT)fc a . 

We have shown that the imposition of Eq. (|6.2|) implies 
the gauge conditions of Eq. (|6.4|l . We now examine the 
reversed question: Does the imposition of the preferred- 
gauge conditions imply that the coordinate description of 
the horizon is the same in the unperturbed and perturbed 
spacetimes? We shall show that the answer is in the 
affirmative. 

Suppose that on the contrary, the perturbed horizon is 
displaced with respect to its unperturbed position. The 
parametric description of the generators is then 

£>, e A ) = z a {v, e A ) + \£ a (v, e A ), (6.6) 

where the vector A£ Q points from a point identified by 
(v,9 A ) on the unperturbed horizon to a point (carrying 
the same intrinsic coordinates) on the perturbed hori- 
zon. We shall show below that if Eq. I|6.4|l holds, then 
£ a must be tangent to the horizon; it can then be de- 
composed as = ak a + a A e a A for some coefficients a 
and a A . If is tangent to the horizon, then it maps a 
point (v,9 A ) on one generator to another point (v',9' A ) 
on another generator. (If 6' A = 9 A the vector links 
two points on the same generator; this happens when 
a A = 0.) Because the mapping preserves the coordi- 
nate labels, this amounts to performing a transformation 
(v',0' A ) — > (v,9a) of the horizon's intrinsic coordinates. 
This transformation can always be undone, and we con- 
clude that £ Q can be made to vanish whenever it is tan- 
gent to the horizon: k a ^ a = => = 0. By showing 
that fc a £ Q = follows from Eq. (|6.4|) we therefore prove 
that Eq. lE3|> implies Eq. (|B~2l . 

According to Eq. I|6.6|l the perturbed basis vectors are 
k a = k a + X^k 13 and e\ = e a A + \£ a e A . The per- 
turbed metric at the new horizon position is g a p{z+\S,) = 
g a p{z + K) + Xh a p(z) = g a/ 3 + X(g a p,-yC + h a p), with all 
fields now evaluated at z, the position of the unperturbed 
horizon. Let now e" stand for any one of the vectors k , 
e A , and e" for the corresponding unperturbed vector. 
The statement g a /3(z + A£)/c Q ef = 0, when expanded in 
powers of A, leads to 

(£^g a + h aS3 )k a e (i a = 0. 

The gauge conditions h a pk a e^ — therefore imply 
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k a e^£^g a p = 0, or 
k a k p £ 



Q 



(6.7) 



These equations come with an immediate interpretation: 
££9a/3 is the change in h a p produced by a gauge trans- 
formation generated by the vector field £ Q ; Eq. (|6.7() in- 
dicates that this transformation must preserve the gauge 
conditions of Eq. i|tj.4fl . 

The first of Eq. (|6.7[1 can be expressed in the form 
(k a t a ) ; pkP - k%k^ a = 0, or 



dv 



(6.8) 



after using Eq. 1)3. 2[) . If we restrict ourselves to situations 
in which the event horizon starts in a stationary state, 
then £ Q = initially and Eq. implies that k a £ a = 
at all times. We therefore have the statement that £ a 
is tangent to the horizon, and the proof that Eq. (|6.4I) 
implies Eq. I|6.2|l . 

The second of Eq. i|6.7|) can be expressed in the form 

|;( e A^) = ~ ( fcQ ^) + 2 ^(fc^ Q ), (6.9) 

where oja was introduced in Eq. fHty . With k a £ a = 
and £ Q = initially, this equation states that e A ^ a = 
at all times. We therefore see that £ a must be directed 
along k a , so that a point (v,8 A ) is necessarily mapped 
to a point (v',6 A ) on the same generator. This simply 
corresponds to a reparameterization v — » v'(v, 9 A ) of the 
generators, and we conclude that £ Q can be set equal to 
zero whenever Eq. 1)6. 4|1 is enforced. 



B. Metric and Lie derivatives 

In the preferred gauge defined by Eqs. H6.2fl - I|6.4[) . the 
perturbed horizon is in the same coordinate position 
as the Kerr horizon, and the perturbed horizon metric 
is Jab = g a pe a A e B = Jab + A/i Q/3 e^e^, having used 
Eq. (|6.3|l . This means that the metric perturbation is 
given by 



Iab 



u a 

n a f3e A e B . 



(6.10) 



The Lie derivative of j\ B in the direction of k a is 
calculated as £^\ B = (ft. a( 3e^e^) ;7 fc 7 . The covari- 
ant derivative of e A in the direction of k a is evaluated 
in Eq. i|2.14[) . and we have £kl A B = h a p- a e\e B k 1 + 
(ujAe B + e A ujB)h a pkP . Using Eq. (|6.4() . we arrive at 



£klAB = Zctab = h aj 3- 1 e A e B k' y . 



(6.11) 



The Lie derivative of ~f AB in the direction of (f> a is 
calculated in a similar way. Here we need the covariant 
derivative of e A in the direction of <fi a , which is evaluated 
as 



in which we expressed <jP as 



refer back to 



Sec. II C — and substituted Eq. I|2.15[) . Gathering the 
results and using Eq. 1)6.10(1 as well as the gauge condi- 
tions of Eq. (|6.4|l . we obtain £$ r y AB = h a (3- 1 e A e B (j) y + 
^ibiA^^c^ ■ We now use the identity c A = cj) D , A — 
T Ac cf> c (which follows from the definition of intrinsic co- 
variant differentiation and the fact that <f> c does not de- 
pend on the horizon coordinates) to manipulate the terms 
involving the horizon connection. We arrive at 

£<i>lAB = K0- n e A e B <jP + cf-ycB + c b1ca- (6-12) 

Recall that the antisymmetric two-tensor cab = —<t>A\B 
was defined and evaluated in Sec. II C; see Eq. 12.1711 . 

We now define a four-dimensional version of this tensor 
with the relation 



c^ = c AB e%e B , 



(6.13) 



where the indices on cab are raised with j AB , the inverse 
of the Kerr horizon metric; this relation is inverted by 
cab = c a pe<\e B . With Eqs. (|6~TU|) and ijETTSj) we have 



C„,i 



c aTcb = c ai h S 0e A e B [y 



(^,CD i j 



Using now the completeness relations of Eq. (|2.12|) and 
the properties c^k 1 = c ai N^ — 0, we obtain c a Jq B = 
cJh 1 pe° A e B and Eq. (|6.12|) becomes 

£<I>1ab = h a j3- 1 e c A e B f < + (cj h li3 +c^ h ja )e A e B . (6.14) 

Finally, the Lie derivative of j AB in the direction of 
t a is calculated as £tl A B — £kJAB ~ ^h£i/>7ab- From 
Eqs. (|6.11|l and l|6.14|l we obtain 

£tlAB = h al 3- 1 e A e B t~< - fi H (c Q 7 ft. 7/ 3 + c^h ia )e A e B . 

(6.15) 



C. Fluxes 

It is a straightforward matter to substitute Eqs. I|6.11[) . 
(|6.14l) . and l|6.15|l into the flux formulae of Eqs. (|4.22|) - 
KZfy . We first obtain 

<J AB £ kl \ B = ±F[k], 
o AB £^i\b = \{F[$\ + G), 
a AB £ tl \ B = i(F[*]-fi H G), 



where 



and 



F[£\ = h af3n k^h^ x e{e Aa c^ A )(e B ^e B ) 



-A;i 



(p AC k a 



kP 



G = /i Q/3;7 fc 7 (c*h x » + c y x h Xli ) {c Aa c A ) (c B ' 3 c B ) . 
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By virtue of Eq. (S!5fr . e Aa e^ = g a ^ + k a N^ + N a k^ , 
and this identity can be used to simplify our expressions 
for F[£] and G. To simplify things further we also write 

and we note that according to Eqs. (|2.13[) and i|2.14[) . 

any tangential derivative of the form £^ 7 £ 7 is necessarily 
proportional to k@. The preceding equation therefore 
becomes 

where p[£] is a proportionality factor that depends on 
the choice of vector £ Q ; for example, p[k] — n and 
p[<p] = ujA<f> A - At this stage it is convenient to supple- 
ment the three gauge conditions of Eqs. Q6.4JI with the 
fourth condition implied by Eq. Ij6.5|) : adopting the radi- 
ation gauge allows us to set h a pkP = on the horizon, 
and therefore to discard all terms of the form fc^/i Q/ 3 ;7 £ 7 . 
This greatly simplifies our expressions for F[£\ and G. 
After carrying out these manipulations we obtain 

F[£] = h a ^h a %i s 

and 

G = 2h a0n k~<c aS h/. 

Substituting these results into Eqs. I|4.22[1 - II4.24[1 wc 

arrive at 

(M) = j,T aP k a t fi dS -0 H j> qdS, (6.16) 

(j) = - j>T al3 k a (j) f:S dS - I qdS, (6.17) 

£{A) = j>T aP k a k^dS, (6.18) 

where 

T a p = -^(h^-ah^.J) (6.19) 

and 

q = -L(h w i a A A "), (6.20) 

167T 

We recall that these results are formulated in the radia- 
tion gauge of Eq. i|6.5[) . And we mention that an alter- 
native expression for q is 

this follows from substituting the (easily-derived) iden- 
tity c al3 = -4> a -,0- K-^N^kp + ^N^ka into Eq. (|6~2U|) 
and simplifying the result. 

The integrals involving T a p in Eqs. I|6.16|l - H6.18|l are 
formally identical to the flux formulae that would be ob- 
tained for a horizon perturbed by a matter field with 



stress-e nerg y tensor T a p (see, for example, Sec. 6.4.2 
of Ref. Hzl, or Sec. 4. 3.6 o f Ref. M)- It is therefore 
tempting to view Eq. (|6.19|l as a definition of an effec- 
tive stress-energy tensor for gravitational radiation cross- 
ing the event horizon. While in general the integrals 
involving q spoil this interpretation, we see that there 
exists an approximate regime in which the interpretation 
is sound: this is the high-frequency regime first inves- 
tigated by Isaacson |3lLl33|. Schematically, T ~ (V/i) 2 
while q ~ KS7h, and the additional derivative ensures that 
T a p dominates over q in the high-frequency limit. And 
indeed, our expression for T a p, as given by Eq. H6.19fl . 
does agree with Isaacson's effective stress-energy tensor. 
It should be noted that the time averaging involved in 
Eq. H6.19fl is different from the spacetime (Brill-Hartle 
[39]) averaging used in Isaacson's construction; but it is 
plausible that the two averaging procedures are recon- 
ciled after T a p is integrated over dS. It should also be 
noted that while Eq. I|6.19l) is formulated in the radia- 
tion gauge of Eq. (|6.5|l . Isaacson has shown that the ex- 
pression is actually gauge invariant in the high-frequency 
limit. 

The flux formulae of Eqs. I|6.16[l - (|6.18|) are not lim- 
ited to the high-frequency regime; they can applied in 
general situations, provided that the metric perturbation 
h a /3 satisfies the gauge conditions h a pk^ — on the hori- 
zon. These formulae could in principle be used in tandem 
with Chrzanowski's metric reconstruction [38l Eol EH E^ 
to calculate the absorption of mass and angular momen- 
tum by a Kerr black hole. But to proceed like this would 
be much more involved than to proceed directly with the 
curvature formalism of Sec. V. The flux formulae could 
also be used in the context of a Schwarzschild black hole, 
but the formulation given here is not optimal and I shall 
refine it in the following section. My main purpose in this 
section was to introduce the preferred gauge (which will 
be used also in Sec. VII) and to establish the preceding 
connection with Isaacson's effective stress-energy tensor 

MM 



VII. METRIC FORMALISM FOR 
SCHWARZSCHILD BLACK HOLES 

In this section I fulfill the promise made in Sec. VI, to 
translate the flux formulae of Eqs. (|4.22() - (|4.24|) into a 
more practical language that involves the metric pertur- 
bations of a Schwarzschild black hole. The key aspects of 
the theory of first-order perturbations of this spacetime 
are summarized the Appendix. 



A. Background spacetime 

The Schwarzschild metric in Eddington-Finkelstein co- 
ordinates {v,r,9,4>) is given by 

ds 2 = -fdv 2 + 2dvdr + r 2 dfl 2 , (7.1) 
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where / = 1 - 2M/r and dfl 2 = n AB d0 A d6 B = d0 2 + 
sin 2 8 d<fr 2 . The subset of coordinates (v,9,4>) is used 
on the horizon; v is a parameter on the null genera- 
tors, and 9 A = (9, </>) are comoving coordinates. In 
the spacetime coordinates (v, r, 9, <f) the basis vectors are 
k a = (1,0,0,0), N a = (0,-1,0,0), e£ = (0,0,1,0), and 
e 4> = (0? 0, 0, 1). The metric of the unperturbed horizon 
is Jab = r+VLAB, where r + = 2M. 



B. Metric perturbation 

The metric perturbation \h a p (denoted 5g a p in the 
Appendix) is cast in the radiation gauge of Eq. I|tj.5[l . 
We therefore impose h a ^k^ = on the perturbed hori- 
zon, which is still located at r — r + — see the discus- 
sion of Sec. VI A. The gauge conditions imply that the 
components h va of the metric perturbations all vanish. 
According to Eq. i|6.10[) the perturbation of the horizon 
metric is 



Iab = h a pe a A e B , 



(7.2) 



where e a A are the background basis vectors. 

The odd-parity sector of the perturbations is described 
by Eqs. (|A.14(1 . (|A.15|) and it involves the functions 
h ™ and h™ of the coordinates (v,r). The combina- 
tions of Eq. i|A.16|) are gauge invariant, and they are 
used in Eq. I|A.18() to form the Regge- Wheeler function 
^Rw( w > r ) 3' which is also gauge invariant. A simple 
calculation shows that near the horizon, 



1 



dhl m 



RW 



2r + dv 



O(f), 



so that ]i2{v,r + ) = 2r + 
parity sector of Eq. I|7.2|) is therefore 



f v ^(v>,r + )dv>. 



The odd- 



sources, K lm = \l(l+l)G lm +0(f) 
\l(l+l)G lm (v,r + ) =JMI + I)r-H% 
parity sector of Eq. i|7.2[l is therefore 



sothatK /m (w,r + ) = 
(v,r + ). The even- 



l.cvcn / 



J2 Z T B (d A )*™(v,r + ), (7.4) 



In 



where Z% 



A ) are the even-parity tensorial harmonics 
introduced in Eq. 1A.7|) . 

The complete perturbation of the horizon metric is 
given by the sum of Eqs. I|7.3j) and 17.4jl . 



Iab = 



£ 

i in 



2X 



AB 



^(v')dv' + Zt B ^(v) 



where we have set f„ v (n') 



¥^(v,r+ 



'ZM 

n AB 



Notice 



Iab 



(7.5) 

w',r+)and¥^(tO = 
that this tensor is tracefree: 
0. The fact that "f AB is related to the integral 
of the Regge- Wheeler function means that this variable 
is rather ill-suited to describe ingoing gravitational radi- 
ation crossing the event horizon; a similar statement is 
made in subsection 5 of the Appendix, about outgoing 
radiation at future null infinity. 

The shear tensor is obtained by differentiating ^J AB 
with respect to v, and 



& AB 



2^*S^(«) + ^S*g&(«) 



(7.6) 



this is also equal to \£kl\ B an( l \ £tl\ B ■ Because the 
spherical harmonics are all proportional to e wn * > , we also 
have 



£-4>1ab 



lm 



2X 



AB 



^(v')dv' 



z ab^zm(v) 



(7.7) 



Iab ( v >° 



l in 



= 2r + > X 



lm j 
AB \ 



¥^(v',r + )dv', 



where X l ™ B 



(7-3) 

) are the odd-parity tensorial harmonics 
introduced in Eq. (|A.8J| . Here and below, the sum over I 
is restricted to Z > 2, and the sum over m extends from 
-I to I. 

The even-parity sector of the metric perturbations is 
described by Eqs. I|A.20(I - (|A.22(I and it involves the func- 
ira j^im ^ an( j Qim Q f tne coordinates (v,r). 



tions h lm 



rr 5 Jr 

The combinations of Eq. (|A.23(I are gauge invariant, and 
they are used in Eq. (|A.24I) to form the Zerilli-Moncrief 
function #£m(^ r ) R 13, 43], which is also gauge invari- 
ant. A simple calculation shows that near the horizon, 



*ZM — 



4r 2 _ 
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1(1 + l)(l 2 + I + 1) dv 
2r + 



K 



Irn 



1 



-l(l + l)G 



lm 



i(i + i 



■K l 



O(f). 



On the other hand, an analysis of the linearized field 
equations near the horizon shows that in the absence of 



C. Fluxes 

It is a straightforward task to substitute the preceding 
results for gab, £klAB, £tjAB, and £^ab into the flux 
formulae of Eqs. (|4.22|) - (|4.24l) . and then to integrate over 
dS = r 2 , sin# d9d<p. The integrations can be carried out 
explicitly with the help of the orthogonality relations of 
Eqs. (|A.11() . i|A.12[l . and we arrive at 



lm 



ZM 



(v)[ 



(7.8) 



and 



(J) 



L^(i_ + i)( l + 2 )(<m) 



64tt 



x U^(v) / ^(v')dv' 



+ *zm(«)*zm(«) 



(7.9) 
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Notice that except for the substitution u — > v, these for- 
mulae are identical to Eqs. I|A.26|) and 1A.27I) . which give 
the rates at which energy and angular momentum are 
transported to future null infinity. Note also that for a 
nonrotating black hole, the first law of black-hole me- 
chanics reduces to (k/8it)(A) — (M). Finally, note that 
although it involves complex quantities, the expression 
for (J) is real; this property follows from the identity 
_ (_2) m ^ r ' m , which is inherited from the spheri- 
cal harmonics, and which is satisfied by both the Regge- 
Wheeler and Zerilli-Moncrief functions. 

Equations (|7.8|l and l|7.9|l give the final form of the 
flux formulae for the case of a Schwarzschild black hole. 
The steps required to compute (M), (J), and (A) are 
therefore these (see also Sec. I D): First, solve the Regge- 
Wheeler and Zerilli Q equations for the functions 
V^hw^' r) and ^zm(v, r) defined in the Appendix, for all 
relevant value of I and m. Second, evaluate the functions 
at r = r + and compute the integral of ^> l ^ v (v,r + ) and 
the derivative of w^ /[ (v,r+). Third, and finally, substi- 
tute these functions into the flux formulae of Eqs. (|7.8|l - 
O and sum over I and m. 

These flux formulae were first presented and used by 
Martel 01 m ms exploration of gravitational-wave pro- 
cesses associated with the motion of a small-mass body 
around a Schwarzschild black hole. Although he arrived 
at the correct results, the derivation of Eqs. I|7.8[l and 
(|7.9|) given by Martel is flawed — it incorporates both 
a conceptual and a computational error. The concep- 
tual error is that Martel based his derivation on Isaac- 
son's effective stress-energy tensor for gravitational waves 
pHl 132] ] , incorrectly assuming that the high-frequency de- 
scription is always applicable near the event horizon of 
a black hole (as it always is near future null infinity). 
This assumption was motivated by the observation that 
for a stationary observer just above the event horizon, 
any incoming gravitational wave would appear highly 
blueshifted. While the observation is of course valid, the 
observer-dependent blueshift does not by itself produce 
a perturbation that satisfies the assumptions underly- 
ing Isaacson's construction — the static Schwarzschild 
coordinates do not form a "steady" coordinate system 
near the horizon. Martel's starting point was therefore 
Eqs. (|6.16(l and . 1 T|> without the integrations over q, 
and this should have led him to the wrong formula for 
the flux of angular momentum (the q integral does not 
contribute to (M), because J7h = for a nonrotating 
black hole). That he nevertheless obtained Eq. (|7.9|l is 
due to a computational error that accidentally compen- 
sated for the absence of the q integral. 



VIII. SMALL-HOLE/SLOW-MOTION 
APPROXIMATION FOR A SCHWARZSCHILD 
BLACK HOLE 

In this and the following section I describe an applica- 
tion of the flux formulae obtained in Sec. V and VII. I 



shall evaluate (M) and (J) in a small- hole/slow-motion 
(SH/SM) approximation in which the ratio M/1Z, where 
M is the black- hole mass and 7Z the radius of curvature of 
the spacetime in which the black hole moves, is assumed 
to be small. I begin in this section with a nonrotating 
black hole, and I will consider the case of a rotating black 
hole in Sec. IX. 

A. The SH/SM approximation 

We imagine a situation in which a black hole of mass 
M is not at rest and isolated, but moves in a spacetime 
that may contain a number of additional bodies. The 
radius of curvature of this external spacetime is denoted 
1Z, and although this may depend on M (if the geometry 
of the external spacetime is significantly influenced by 
the black hole), we assume that 

M/1Z <C 1 (small- hole/slow- motion approximation). 

(8.1) 

More precisely, we assume that M is much smaller than 
all of 7Z, C, and T, where L is the scale of inhomogeneity 
in the external universe, and T is the time scale over 
which changes occur in the external universe. To simplify 
the notation we take 1Z, C, and T to be of the same 
order of magnitude. (These quantities, and many of the 
concepts used throughout Sees. VIII and IX, are defined 
precisely in Thorne and Hartle [33j ; the reader is referred 
to this paper for details.) 

Near the black hole the spacetime resembles closely 
the spacetime of an isolated black hole: the gravitational 
field is strongly dominated by the hole's contribution, 
and the influence of the external universe is weak. But 
the hole is not truly isolated, and it is slightly distorted 
by the tidal gravitational field supplied by the external 
universe. As a result of this interaction, the hole's mass 
and angular momentum change with time, and we wish 
here to calculate these changes. 

When viewed on the large scale 1Z, the black hole occu- 
pies a very small region of the actual spacetime, and this 
region can be idealized as a world line 7 in the external 
spacetime. Let u a be the (normalized) tangent vector to 
this world line, and call this the four- velocity of the black 
hole in the external spacetime. It can be shown that to 
a very good degree of accuracy, the motion of the black 
hole is geodesic in this spacetime [T9L El liH Eol EtI l48| . 
Let e" (with the index a running from 1 to 3) be a set of 
orthonormal vectors attached to 7; let these vectors be 
orthogonal to u a and choose them to be parallel trans- 
ported on the world line. The tetrad (u a ,e") defines a 
reference frame in a neighborhood of 7, and we shall call 
this frame the local asymptotic rest frame of the black 
hole in the external spacetime. 

We assume that the Ricci tensor of the external space- 
time vanishes in a neighborhood of 7, so that no mat- 
ter will appear in the vicinity of the black hole. The 
curvature of the external spacetime in this neighbor- 
hood is therefore described entirely by the Weyl tensor 



22 



C a -yps- The Weyl tensor evaluated on the world line can 
be decomposed in the tetrad (w Q ,e"); we write, for ex- 
ample, C 0c0d (v) = C ai ps{n)u a eju p e & d and C Qcbd {v) = 
C a -)ps(l)u a e2e^e d , with v denoting proper time on 7 — 
it will later be identified with an advanced-time coordi- 
nate on the black-hole horizon. It is easy to show that 
the frame tensors j3^| 

£ab(v) = C aOb(v), B ab {v) = \e f C cdW {v) , (8.2) 

where e a ;, c is the permutation symbol (all frame indices 
are lowered and raised with 8 a \y and its inverse, respec- 
tively), are symmetric and tracefree, and that their com- 
ponents comprise all ten independent components of the 
Weyl tensor. These frame tensors are the tidal gravita- 
tional fields supplied by the external universe, and these 
are responsible for the tidal distortion of the black hole. 

It will be convenient to promote the tidal fields £ a & and 
B a b to four-dimensional spacetime tensors. We therefore 
define 



£a(3 — £ab£a e % 



B n 



^ab e a e fii 



(8.3) 



where e a a = 6 ' g a ^e^. It is not difficult to show that 
these tensors are also given by 



and 



Bc.fi = \u»e^ s C, 



(8.4) 



(8.5) 



where the Levi-Civita tensor e^ au p and the Weyl tensor 
evaluated on the world line 7. 
As an example of a SH/SM situation, consider a black 
hole of mass M on a circular orbit of radius b in the grav- 
itational field of an external body of mass M ext . The 
radius of curvature of the external spacetime at the posi- 
tion of the black hole is such that TZ~ 2 ~ (M + M ext ) /b 3 , 
and we have 



71/ 

n 



M 



M + M P 



-V 3 



V = 



M + M cxt 



(8.6) 



where V is a measure of the hole's orbital velocity. There 
are many ways by which M/TZ can be made small. One 
way is to let M/M cxt < 1; then M/TZ will be small 
irrespective of the magnitude of V. This is the small-hole 
approximation, which allows the small black hole to move 
at relativistic speeds in the strong gravitational field of 
the external body. Another way is to let V <C 1; then 
M/TZ will be small for all mass ratios. This is the slow- 
motion approximation, which allows the slowly-moving 
black hole to have a mass comparable to (or even much 
larger than) M ex t ■ These two limiting approximations are 
special cases of the fundamental requirement that M/TZ 
be small; we therefore call the approximation M/TZ <C 1 
the SH/SM approximation. 



B. Metric of a tidally distorted black hole 

My considerations thus far have been general, and they 
apply to rotating as well as nonrotating black holes. I 
now specialize to nonrotating black holes. 

The metric of a Schwarzschild black hole immersed in 
an external universe can be obtained by solving the Ein- 
stein field equations. Because the tidal potentials scale 
as (r/TZ) 2 <C 1, where r is a measure of distance from the 
black hole, it is sufficient to linearize the equations with 
respect to the Schwarzschild solution; this is a standard 
application of black-hole perturbation theory. Explicit 
forms for the metric were obtaine d by Manasse 0] , Alvi 
[49|, Detweiler |50|, and Poisson [5l|, and I summarize 
their results here. I follow the description of Ref. [5l| . 
but I switch from the retarded coordinates (u, r, 6 A ) used 
there to a set of advanced coordinates (v, r, 9 A ) which are 
well behaved on the event horizon; the expressions for the 
perturbed metric are identical, except for the correspon- 
dence du — > — dv. 

The metric takes the form of an expansion in powers 
of r/TZ, but it is correct to all orders in M/r. It is given 
by 



9vA 
9AB 



-f{l + r z f£*) + 0{r i /TZ A ) 



-yf(£* A +B* A ) + o(ryn% 



(8.7) 
(8.8) 

(8.9) 



r 2 n AB V 

+ 0(r /TZ 3 ), 



2M 2 . , 
1 s- \£ 



AB "T" &AB 



(8.10) 



where / = 1 — 2M/r. The irreducible tidal fields are 
defined by 



£* = 5> m r m , 

m 

m 

m 

m 



A J 



A ' 



B 



AB 



ABi 



(8.11) 
(8.12) 
(8.13) 
(8.14) 
(8.15) 



where Y m , Y. 



A- 



^ AB ' A A' 



and X AB are the real spher- 
ical harmonics of degree I = 2 that are introduced in 
subsection 2 of the Appendix, and 



(8.16) 
(8.17) 
(8.18) 

(8.19) 
(8.20) 



£0 


— £33 — 


-(fn + fiaa), 


Sic 


= 2£i3, 




£is 


= 2£ 23 , 
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- \^ 


- £22), 


£is 


— £12, 
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with corresponding relations defining B m . (Notice that a 
typographical error contained in Ref. |5l| is hereby cor- 
rected.) 

In the limit M/r — > (keeping r/TZ fixed), the metric 
of Eqs. H8.7[) - H8.10[) becomes the metric of the external 
spacetime expressed as an expansion in powers of r/TZ 
about the timelike geodesic 7. In this limit the inter- 
pretation of v as proper time on the world line becomes 
precise, and £ a b(v), B a b(v) are recognized as frame com- 
ponents of the Weyl tensor evaluated on 7. In the limit 
r/TZ — ► (keeping M/r fixed), the metric of Eqs. (|8.7fl - 
(|8.1UI) becomes the metric of an isolated Schwarzschild 
black hole expressed in ingoing Eddington-Finkelstein co- 
ordinates; there is no notion of a world line 7 in this 
limit. For small values of r/TZ and arbitrary values of 
M/r, the metric of Eqs. H8.7fl - H8.1l)fl describes a black 
hole distorted by the tidal gravitational fields supplied 
by the external universe. 

The metric perturbation h a p defined by Eqs. (|8.7[) - 
(|8.1()|l satisfies the conditions h a f3k a k f3 = h a/ 3k a e^ = 
at r = r + = 2M, because h vv and h vA are all propor- 
tional to / = 1 — 2M/r. The metric perturbation there- 
fore satisfies the preferred gauge conditions of Eq. 16.41) . 
It can also be checked directly from the metric that the 
hypersurface r = 2M is null, and that its generators move 
with constant values of 9 A . The parameter on the gener- 
ators is v, and a short calculation reveals that the surface 
gravity k is equal to its Schwarzschild value (4M) _1 up 
to a fractional correction of order (M/TZ) 3 . 



C. Odd-parity contribution to shear 

Although the metric of the tidally distorted black hole 
is already expressed in the preferred gauge, it is safer 
(and as it turns out, necessary) to calculate the shear 
tensor <Jab by first obtaining the gauge-invariant Regge- 
Wheeler [3 and Zerilli-Moncrief 0, \l3L I43I ] functions; the 
relation between these quantities is given by Eq. I|7.6[) . 
We begin here with the odd-parity piece of the shear 
tensor. A description of this sector of the metric pertur- 
bations is provided in subsection 3 of the Appendix. 

When / = 2 the odd-parity perturbations can be ex- 
panded as 



and 



h iA = Y,K{v,r)X<Z(9 A ) 



h AB = Y. h 2( V > r ) X AB(0 J 



The combinations 



1,™ 1 



.21) 



.22) 



(8.23) 



are gauge invariant, and the Regge- Wheeler function is 
defined by 

1 



*Kw = ~r' l h?. (8.24) 



Com parison of E qs. ED, gjUjl with Eqs. (pO|l . 
ijOty using Eqs. fsHJl . ^J^i reveals that h™ = 

\r A B m . Equation 



1^3 



B m , up to 



-~\r 3 fB m , h? = 0, and h" 2 
then gives h™ — — ^r 3 fB m and h™ 
smaller terms proportional to dB m /dv ~ TZ~ 3 . From 
Eq. (|8.24|) we obtain = 0. This curious result leads 
to the conclusion that the metric of Eqs. Ij8.7|l - i|8.10|l is 
not sufficiently accurate to calculate the Regge- Wheeler 
function, and therefore the shear tensor. 

Fortunately, the Regge- Wheeler equation i|A.19|l is suf- 
ficiently simple that it can be solved directly. Assuming 
(as we shall verify below) that derivatives of ^rw with re- 
spect to v can be neglected compared with spatial deriva- 
tives, the Regge- Wheeler equation for I = 2 reduces to 



r(r-2M)^+2M^-- 
dr A dr 



M 

■61 1 

r 



^ RW = 0. (8.25) 



The solution that is well behaved at the horizon is vI/rw <x 
r 3 , and to produce the correct metric perturbation we 



write 



(8.26) 



where the overdot indicates differentiation with respect 
to v. While B m (v) scales as TZ~ 2 , its time derivative scales 
as TZ~ 3 and w-derivatives of the Regge- Wheeler function 
are indeed much smaller than its spatial derivatives. 

To see that Eq. I|8.26l) is indeed the correct solu- 
tion to the Regge- Wheeler equation, we reconstruct the 
metric perturbation in the preferred gauge and show 
that it agrees with the odd-parity sector of Eqs. (|8.7|l - 
(18.101) . We note first that according to Eq. I|8.24|l . 
^rw — r {hv + fh>r)> where we have removed the la- 
bel m for simplicity. On the other hand, the field equa- 
tion QA. 17(1 implies d v h r + d r (h v + fh r ) = 0. Solving 



these equations yields h v 



1^.3 
3 



r A fB and h r 



1„3 



r 3 B, 



up to smaller terms involving B. The actual metric per- 
turbations are then recovered by using Eq. I|8.23|l along 
with the gauge condition h r = 0. The equation for h r 
gives {\d r — f _1 )/i2 = hr, and this differential equa- 
tion has /12 = \ rA B as solution; this agrees with our 
previous expression. Finally, the equation for h v gives 



hd v h 2 — h v , and we obtain h v 



— _ 1 „3 



r 3 fB up to a 



smaller term involving B; this also agrees with our pre- 
vious expression. We conclude that the Regge- Wheeler 
function of Eq. (|8.26|l is indeed compatible with the met- 
ric of Eqs. (fH~7|> — (fHTTUf) . 

Substituting Eq. {0| into Eq. (JOl) we obtain 

<b(v,0 A ) = ~riJ2^ X AB(0 A )- (8.27) 



According to Eq. I|8.15|l . this can also be written as 

AB — 12' + a AB- 
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D. Even- parity contribution to shear 

We turn next to the even-parity sector of the metric 
perturbations; the reader is referred to the description 
given in subsection 4 of the Appendix. 

When I = 2 the even-parity perturbations can be ex- 
panded as 

hij = ]T^>,r)y m (0, (8.28) 

m 

h iA = £ir>,r-raO, (8-29) 

m 

h AB = r 2 Y / G r "(v,r)Z% B (e A ), (8.30) 

m 

where we have incorporated our knowledge that K m = 
3G m ; this follows because the metric perturbation of 
Eq. (|8.1U|) is tracefree [refer also to the discussion that 
precedes Eq. (|7.4|) ]. The combinations 



h m = h m 



2e 



K m = 3G m - -r x e^, 



where e™ — j™ — ^r 2 G™,a,re gauge invariant, and the 
Zerilli-Moncrief function [j, \13L |43 | is defined by 



v& m = _ 

*zm — o 



A v 13 



.32) 



where A = 4 + 6M jr. 

Com parison of E qs. ifOfy-ljOty with Eqs. l|8~7)l 
ifOty using Eqs. |HHIJ— (JHUBJ reveals that h™ v = 
-r 2 f 2 £ m , = -y/£ m , G m - -±r 2 (l - 2M 2 /r 2 )£ m , 
and if m = 3G m , with all other components vanishing. 
From this information the Zerilli-Moncrief function can 
be computed straightforwardly (no need to solve the Zer- 
illi equation directly), and its value on the horizon is 
found to be 



*2 M (v,r + ) = -~r 3 + £ m (v). 
Substituting this into Eq. I|7.6[) we obtain 



.33) 



.34) 



According to Eq. Q8.13p . this can also be written as 



AB 



~ 12 ' +°AB- 



E. Shear and Weyl tensor 

The sum of Eqs. I|8.27|l and l|8.34|l gives the complete 
shear tensor, 

ctab{v,9 a ) = -^rlJ2[^( v ) Z A B (0 A ) 

m 

+ B m (v)X% B (0 A j\, (8.35) 



and we observe that this 2-tensor is properly tracefree. 

We can use Eqs. I|8.35|l and Eq. (|4.7|) to calculate the 
Weyl tensor Cab on the horizon. We may neglect time 
derivatives and write Cab — no ab , which gives 



Cab = 



£ m Z 



m^AB 



B m X 



AB 



(8.36) 



This expression shows that the (dimensionless) Weyl cur- 
vature on the horizon is of order (M/1Z) 3 , which is a fac- 
tor M/1Z <C 1 smaller than the asymptotic value of the 
Weyl curvature (for r 3> M). 

For future reference we calculate the Weyl scalar \f 
from Cab] this is defined by Eq. I|5.6[l and related to the 
Weyl tensor in Eq. 1)5. 5fl . From this equation and the 
properties of the vectors e A we infer that ^> — Cab& A & B - 



With Eqs. (|A~9jl and l|A~T0)l we can relate the tensorial 
harmonics Z™ B and X™ B to the spin-weighted spherical 
harmonics ±2Y m . (Please note that the vectors e A used 
in the Appendix are rescaled versions of the vectors used 
here: e A = e A /r + ; in the notation of the Appendix we 



.31) have ^ = Cab£ e / r + ) Simple algebra then gives 



V{v,6 A ) 



V6 
24 ? 



(v)-iB m (v) 2 Y m (9 A ). 



(8.37) 

This reveals that \& = 0(M/1Z 3 ) on the horizon, while 
* = 0(1/K 2 ) asymptotically (for r > M). This shows 
once more that the Weyl curvature on the horizon is sup- 
pressed by a factor M/1Z <§; 1 with respect to its asymp- 
totic value. 



F. Fluxes 

The shear tensor of Eq. 1)8. 35|) is equal to \£u1ab — 
\£tl\ B , an d the metric perturbation r )\ B {v 1 9 A ) can be 
obtained by direct integration with respect to v. Differ- 
entiation with respect to then gives £<f>JAB> this can 
be worked out by using the explicit forms for Z™ B and 
X™ B gathered from subsection 2 of the Appendix. Fi- 
nally, these results can be substituted into the flux for- 
mulae of Eqs. (|4.22fl - l|4.24[l . and integration over dS — 
r 2 , sin 9 d9d(f> is readily carried out using the known an- 
gular dependence contained in the spherical harmonics. 
A straightforward computation yields 



(M) = 8 -^(3£ 2 +£l + £l+A£l + A£l 
+ 3B 2 +B 2 c + B 2 ls +4B 2 2c + 4B 2 s ) 



and 



(j) 



8M e 
^t5~ 



(£ic£is — £\s£\c + 8£2c£2s — 8£2 S £2c 



B lc B ls - B ls B lc + 8B 2c B2s - 86 2s i3 2c ), 



where £ m and B m are the harmonic components of the 
tidal gravitational fields introduced in Eqs. I|8.16[) - I|8.2(J[I . 
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These results can be expressed in terms of invariants 
formed from £ ao and B a b, the components of the tidal 
fields in the local asymptotic rest frame of the moving 
black hole. We also need the derivatives of these fields 
with respect to v (denoted with an overdot), and the 
unit vector s a = (0, 0, 1) that points in the direction of 
the third coordinate axis. (This direction is preferred 
because the angles 6 and <f> refer to it.) In terms of these 
quantities, the previous expressions become 



( M) = 16M6 ^ *•><■ 



45 



(8 ab £ ab + B ab B ab ) 



and 



(j) 



32M 6 
45 



e acd (£ a b £ bc + B\B hc ) 



bc\ s d 



(8.38) 



(8.39) 



where e ac d is the three-dimensional permutation symbol. 
We also have (k/8tt)(A) = (M). We note that when ( j) 
is expressed in the covariant form of Eq. (|8.39|l , what is 
actually meant by (J) is the rate of change of the compo- 
nent of the angular-momentum vector in the direction of 
s a ; in three-dimensional vectorial language appropriate 
in the local asymptotic rest frame, (J) = (J a )s a , where 
J a is the vectorial angular momentum. 

Alternatively, the flux formulae can be expressed in 
terms of the spacetime tensors of Eq. (|8.4ll and 18.511 . 
The translation is effected by Eq. 1)8. 3fl and the identity 

£abc — ufl £fj,apyea e b e c> where £ Ma /3 7 is the Levi-Civita 
tensor. We find 



1 fiM 6 

(M) = -^(£ a0 £«P + B a/3 B af} ) (8.40) 



and 



(j) 



32Af f 
45 



u»e miS (£y^ + B a ^)s 5 , (8.41) 



where s a = s a e" is a unit spatial vector, and £ a p = 
£af}-,^u tl , B a p = B a f3-^u^ are the proper-time derivative 
of the tidal gravitational fields. We recall that in this 
SH/SM description, all vectors and tensors refer to the 
spacetime of the external universe in which the black 
hole moves. From Eqs. H8.38|) - l|8.41|) we gather that (M) 
scales as M 6 /1Z 6 , while (j) s c ales as M 6 /K 5 . To the best 
of my knowledge, Eqs. H8.38[) ~ H8.41|) have never appeared 
before in the literature. 



G. Comparison with Thorne, Hartle, and Zhang 

The rate of change of angular momentum for a general 
body interacting with a tidal gravitational field was cal- 
culated, in the regime M/1Z <C 1, by Thorne and Hartle 
|33|: they obtained the expression 



be 



dc 



-,J b d ® 



dc 



.42) 



where M a b is the body's mass quadrupole moment, while 
J a b is its current quadrupole moment (both defined in 
terms of the structure of the gravitational field outside 
the arbitrary body). Zhang [53], on the other hand, cal- 
culated the rate at which the body changes its mass; he 
obtained an expression equivalent to 



(M) = -(M ab £ 



?ab 



-JabB 



bub 



(8.43) 



We wish to show that our previous results are compatible 
with these expressions. 

An isolated Schwarzschild black hole is spherically 
symmetric, and its intrinsic quadrupole moments vanish: 
M a b = Jab = 0. But a black hole immersed in an external 
universe is tidally distorted and therefore acquires nonva- 
nishing moments. It is easy to see that Eqs. (|8.38|) . (|8.39|) 
are compatible with the general results of Eqs. H8.42fl . 
(|8.43l) if the tidally-induced quadrupole moments of a 
nonrotating black hole are given by 



_32Af 6 - 

Mab — 77 — O a b 

45 



and 



Jab 



8A/ £ 



Bab 



(8.44) 



.45) 



These scale as M 3 (M/1Z) 3 , and they both involve the 
rates of change of the tidal gravitational fields. This 
is a rather surprising result, as one would expect the 
quadrupole deformation of a tidally-distorted body to 
be proportional to the tidal gravitational field itself, in- 
stead of its time derivative. But the time derivative is 
present, and its origin can be traced back to Eq. H8.37fl : 
the Weyl curvature at the horizon is proportional to the 
time derivative of the asymptotic curvature. Since it 
is the horizon curvature that produces the black-hole 
distortion, this explains why a time derivative enters 
Eqs. (|Q4) and j05)l . 



H. Black hole in a circular binary: Slow-motion 
approximation 

If we specialize to a slow-motion situation, the tidal 
gravitational fields of the external universe can be ap- 
proximated by 



Sab — &,ab 



Bab ^ 0, 



where $ is a Newtonian potential. For concreteness we 
take the Newtonian field to be produced by an external 
body of mass M ext located at r ext (i) relative to the sys- 
tem's center of mass. Then &(x) = — M cxt /\x — r cxt |, and 
we exclude the contribution — M/\x — r \ from the black 
hole because this does not produce a tidal field at the 
position r(t) of the black hole. Also for concreteness we 
take the orbit to be circular, and we let b = \r — r ext | be 
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the constant relative separation between the two bodies. 
The orbital angular velocity is 



n = 



M + M e: 



b 3 



.47) 



and the relative position vector is r — r cx t = P — 
&(cos ttt, sin fit, 0) = bp(t). The relative velocity vector 
is V = 6f2(— sinfii, cosftt, 0) = bflcj)(t). For simplicity 
we align the spin vector in the direction of the orbital 
angular momentum: s = (0,0, 1). 

Using this information we calculate £ a b(t) — 
(M ext /b 3 )(S ab - 3p a Pb) and 



3M cxt f> . 



b 3 



\Pa<t>b + 4>aPb) 



Substituting this into Eq 
32 



gives 



5 V M + M cxt 



M 



V 



18 



(8.48) 



where rj = MM cxt /(M + M ex t) 2 is a dimensionless 
reduced-mass parameter and 



V 



M + Mr, 



< 1 



(8.49) 



is the relative orbital velocity. The rate of change of the 
hole's angular momentum can be obtained directly from 
this and the rigid-rotation relation (J) = f2 _1 (M); this 
gives 

These results agree (in a limit of no black-hole rotation) 
with earlier expressions obtained by Alvi |18| . In the 
regime M cxt <C M they also agree with earlier results 
derived by Poisson and Sasaki 



Black hole in a circular binary: 
approximation 



Small-hole 



We now allow the black hole to move rapidly in the 
strong gravitational field of another Schwarzschild hole 
of mass M oxt ; to comply with the SH/SM condition 
M/TZ <C 1 we now impose M/M ext <C 1, as was discussed 
in Sec. VIII A. Once more we choose the orbit to be circu- 
lar. In the standard Schwarzschild coordinates (t, r, 0, 4>) 
used in the background spacetime of the large black hole, 
the orbital radius is b and the four-velocity of the small 
hole is u a = 7(1,0,0,0), where 7 = (1 - 3M cxt /6)" 1/2 is 
a normalization factor and 



b 3 



.51) 



is the angular velocity. We again align the spin vector 
in the direction of the orbital angular momentum, so 



that s a = (0,0,-1/6,0). Calculation of £ a p, B a p us- 
ing Eqs. (|8.4|) . and substitution into Eqs. I|8.40|) . 
(|8.41|) gives 



(M) = - 



32 ( M \\ rl8 (l-V 2 )(l-2V 2 ) 



and 



(j) = 



5 \M e: 



32/ M 



V 



5 \M e: 



M BXt y 



15 (1-V 2 )(1-2V 2 ) 



(1-3V 2 ) 2 



.52) 



(8.53) 



where V = Af oxt jb < is a measure of the hole's 

orbital velocity. These results agree with those of the pre- 
ceding subsection in a common domain of validity defined 
by M <C M ext and V< 1. To the best of my knowledge, 
the results of Eqs. (|8.52|l and l|8.53|l . complete with all- 
order relativistic corrections, have never appeared before 
in the literature. 



IX. SMALL-HOLE/SLOW-MOTION 
APPROXIMATION FOR A KERR BLACK HOLE 

In this section I apply the SH/SM approximation in- 
troduced in Sec. V III A to the flux formulae derived in 
Sec. V D, Eqs. H5.23[H5.29jl . I wm proceed much as in 
Sec. VIII, except that I will deal with curvature pertur- 
bations — and the Teukolsky equation — instead of 
metric perturbations. 



A. Flux formulae in the SH/SM approximation 

We begin by isolating, in Eqs. (|5.27|l and (|5.28l) . the 

terms for which m = 0: 



(M) 



1 



(r 2 + +a 2 ) / (|$°| 2 )sin0d0 



ntt H J ($™$ m - $™$ m ) sin e d9 



(9.1) 



(j) 



4k 



_ <f>™$™) sin 6 d6. (9.2) 
We recall the definitions 

roc 

= e KV / e - (K - mnH >'y m (v',6)dv', (9.3) 



where 



e mllH " V m (v',6)dv', 



imtp 



(9.4) 



(9.5) 
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is the horizon Weyl scalar introduced in Eq. I|5.6[l . 

To see how we may specialize Eqs. and (19 .4|) to 

the SH/SM approximation, consider first an integral of 
the form F_(v) = J^_ oo e luJV f(v')dv', and suppose that 
f(v') varies on a time scale r that is large compared 
with (We also suppose that / switches off suffi- 

ciently rapidly in the infinite past so that the integral 
converges.) Then F_(v) can be evaluated by successive 
integration by parts, each integration generating a rel- 
ative correction of order e = (wt) _1 <C f. To leading 
order, F-(v) = -iuj- 1 f(v)e^ v [1 + 0(ie)}. Consider next 
an integral of the form F + (v) = f e~ Au /(V) dv' , where 
the real part of A is assumed to be positive; here we sup- 
pose that e' = (At) -1 <C 1. Integration by parts in this 
case leads to F + (v) = X' 1 f(v)e- Xv [l+0(e')]. These sim- 
ple manipulations allow us, within the stated conditions, 
to approximate the integrals by local expressions. This 
is the technique we shall employ to evaluate Eqs. (|9.3|) 
and 

In this way we obtain 



$T = 



* m (w,6»)e imnH1 



K — imD, 



H 



1 + 



1 



(k — imfln)T 



and 



1 + 



1 



where r is the time scale associated with changes in 
fy m (v, 0). To see how the conditions kt > 1 and Oht ^> 
1 relate to the SH/SM approximation, we first recall that 
changes in \t(u, r +1 9, ip) are governed by processes taking 
place in the external universe, so that t ~ 1Z. We also ex- 
press k and Oh in terms of the black-hole mass M and its 
dimensionless rotational parameter \ = a/M = J/M 2 : 



X 



2M(1 + Vl-X 2 ) 



H = 



X 



2M(1 + y/T^j?) 



(9.6) 

we recall that x is limited to the interval < \ < 1 and 
that r + = M(l + y/l — x 2 )- In orders of magnitude we 
have k ~ 1/M and Oh ~ x/M, and to achieve kt 3> 
1 and H r > 1 we need M/K < 1 and M/ft < x, 
respectively. The stronger condition is 



M/n « x, 



(9.7) 



and we take this to be the precise statement of the small- 
hole/slow- motion condition when we deal with rotating 
black holes. Notice that by virtue of Eq. 1)9. 7|) . the no- 
rotation limit x —> will be inaccessible in our analysis; 
this case was treated separately in Sec. VIII. We shall 
write our previous results as 



k — imO 



H 



1 + 0(M/K) 



1 + 0{M/K) 



(9.8) 
(9.9) 



with the understanding that the error terms are really of 
order Mj (x&) , and therefore small by virtue of Eq. (|9.7I) • 
For the remainder of this section we assume that x is of 
order unity, and we allow ourselves to lose sight of this 
distinction. 

Substituting Eqs. I|9.8|l and 1(9. 9|) into Eq. H9.ll) reveals 
that each m/0 term vanishes to leading order in M/K; 
what remains is 



(M) 



2k 2 



(\V°(v, 6)\ 2 ) sm8d8 + 0(M 5 /K 5 



(9.10) 

The scaling of the error term follows from the facts that 
each contribution to a m / term is of order (M/K) 4 , 
but that the cancellation suppresses this by a factor of (at 
least) M/K. An a priori estimate of the surviving term 
in Eq. ||<HU|) indicates that it is of order (M/K) 4 , but 
we shall see that it is in fact of order (M /K) 6 . Inserting 
Eqs. §D| and iJOJ into Eq. gives 



(J) = - 



20 T 



a 2 x , 



(« 2 



TO 2 ^)- 1 



x / (\^> m (v,8)\ 2 )sm8d8, (9.11) 



and this is of order M 5 /K 4 



B. Weyl scalar: asymptotic values 

To proceed further we must compute the functions 
^ m (v,8) that enter into the simplified flux formulae of 
Eqs. (|9.1U|) and l|9.11|l . This requires solving the Teukol- 
sky equation for the Weyl scalar ipo(v,r,8,ip), with ap- 
propriate boundary data provided by the conditions in 
the external universe. In this subsection I specify these 
boundary conditions; in the next subsection (Sec. IX C) 
I tackle the integration of the Tcukolsky equation and 
construct ty m {v, 6) from the solution. My presentation 
in these two subsections will not stray very far from what 
is contained in Sec. 3.5 of the book by D'Eath |19|: and 
my end results will be equivalent to his. 

The function ipo(v, r, 9, ip) we shall work with is 



-Cl 1/3S k a m-«k l3 m s 



(9.12) 



where ^ is the perturbation of the Weyl tensor, while 
k a = k a (K) and m a = m a (K) are members of Kinners- 
ley's null tetrad 0,113 ■ The relation between ip = ipo(K) 
and *£(v, r+, 8, ip) is given by Eq. 1)5.6(1 . 

We wish to calculate i/)q in the asymptotic regime 
r r + , assuming that r is still much smaller than 
1Z, the radius of curvature of the external spacetime. 
The asymptotic values will be constructed from £ a b(v) 
and B a b(v), the tidal gravitational fields introduced in 
Sec. VIII A — Eq. ((8.2(1 : recall that lower-case Latin 
indices refer to the hole's local asymptotic rest frame, 
and that the hole's angular-momentum vector is directed 
along the third coordinate axis. 
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In the asymptotic regime r r + the coordinates 
(v,r,9,ip) are easily related to a set of Lorentzian coor- 
dinates (t, x, y, z) that are adapted to the frame (u a , e"); 
the relations are t — v — r, x = r sin 9 costp, y — 
r sin#sin-0, and z — rcos8. In this regime the null vec- 
tor k a can be decomposed as k a ~ u a + r a , where u a is 
the hole's velocity vector in the external spacetime, and 
r a is a spacelike vector that points radially outward. In 
the asymptotic Lorentzian coordinates (t, x, y, z) we have 
u a = (1,0,0,0) and r a — (0, sin 9 cos ip, sin 9 simp, cos 9). 
In the limit we also have m a ~ 2~ 1 / 2 (0,cos6 l cos?/' — 
i sini/;, cos 9 sin -0 + i cos^, — sin 9). 

We have seen in Sec. VIII A that in the vicinity of the 
black hole (the region r<7?, which includes the asymp- 
totic region r ^> r+), the Weyl tensor of the external 
spacetime can be decomposed into the symmetric, trace- 
free fields £ a b and B a b- In the asymptotic coordinates 
(t, x, y, z) the decomposition is given by Ctatb = £ab, 

Ctabc — —£bcd,B d a i an d Cacbd = $ab£cd + 5 c d£ a b ~ ^ad^bc ~ 

5b c £ad- In these relations the Weyl tensor, and the tidal 
gravitational fields £ a b and B a b, are evaluated on the 
black hole's world line in the external spacetime; they 
are functions of t (or v) only. 

According to Eq. I|9.12|) , the asymptotic value of tpo is 
— C ai f3s(u a + r a )rn 1 (u^ + r /3 )m 5 , with C ai /3s denoting 
the Weyl tensor of the external spacetime evaluated on 
the hole's world line. Using the information provided 
in the preceding two paragraphs, we obtain the explicit 
expression 



-2£ ab m a m b + 2r a m b e abc B c 



(9.13) 



where r a and m a are the spatial components of the vec- 
tors r a and m a , respectively. The angular dependence 
contained in these vectors is encoded in spin-weighted 
spherical harmonics of degree / = 2 (see subsection 2 of 
the Appendix for a definition). It is convenient to intro- 
duce a set given by 



2 V 2 °(M) = --sin 2 9, 
2 Y 2 ±1 (9,^) = -sin6'(cos(9Tl)e ±iV ', 



2 Y 2 ±2 {9^) = 



-(l=p2cos( 



cos 2 9)e ±2 ^. 



(9.14) 
(9.15) 
(9.16) 



then it is straightforward to show that Eq. (|9.13|l is equiv- 
alent to 



Wo 



(9.23) 



This is the asymptotic value of the Weyl scalar 
ipo{v,r,6,ij)) in the regime r+ < r < K, expressed in 
terms of the tidal gravitational fields £ a b{ v ) and B a b(v). 



C. Teukolsky equation 

To relate \D ,m («, 0) to the asymptotic value of tpo ob- 
tained in Eq. (I9.23|l it is necessary to solve the Teukol- 
sky equation [2] for s — 2 and 1 = 2. Because the v- 
dependence of the solution enters through the tidal grav- 
itational fields £ a b(v) and B a b{v), and because this depen- 
dence is slow (time scale of order 1Z), it is actually suffi- 
cient to integrate the time-independent Teukolsky equa- 
tion. We therefore write 

Mv, r,M) = -£ [%(«) + iPm(v)] R m (rhY 2 m (9, V), 

ra 

(9.24) 

where R m (r) is a radial function normalized so that 
Rm(r ^> r + ) ~ 1; this function must be a solution to 
Eq. (2.10) of Teukolsky and Press 0, in which we set 
u) = 0. 

The explicit form of the radial equation is 



x(l 



where 



dx 2 



3(2a; + l) + 2im7 
2a; + 1 



d 

dx 



+4im~f 



x(l + x) 



r 



x = 

r + — r_ 

is a new independent variable, and 

a 



7 



r 



+ 



This set is not normalized; we have instead / 1 2 ^2° | 2 dtt = we recall that r± = M ± V 1 M 2 



R m {x) = 0, (9.25) 
(9.26) 
(9.27) 

The relevant solu- 



24tt/5, / ^Y^dVL = 167r/5, and J \ 2 Y 2 ±2 \ 2 dVl = 4tt/5. 
If we also introduce 



and 



olq = £11 + £22, 
a±i = £13 T i£23, 
a±2 = £11 - £22 T 2i£i2 



/3 = B u + B a2 , 
Ati = B 13 TiB 2 3, 
P±2 = B11 - B22 T 2iBi2, 



(9.17) 
(9.18) 
(9.19) 



(9.20) 
(9.21) 
(9.22) 



tion to Eq. I|9.25|l is 

R m (r) = A m x- 2 {1 + x)- 2 F(-4, 1; -1 + 2im Tl -x), 

(9.28) 

in which the hypergeometric function is actually an ordi- 
nary polynomial of degree 4 in the variable — x. Equation 
(|9.28|) is essentially Eq. (5) from Ref. , and the super- 
ficial difference is attributed to the fact that Alvi works in 
Boyer-Lindquist coordinates instead of our Kerr coordi- 
nates. This is also Eq. (3.7) in Chapter VI of Teukolsky's 
Ph.D. dissertation |53|, and Eq. I|9.28|l is equivalent to 
Eq. (3.5.7) of Ref. [lJj- The constant A m must be cho- 
sen so that the radial function approaches unity when 
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x — > oo; a simple calculation shows that it must be given 
by 

A m = --rnyil + imj)(l + 4to 2 7 2 ) (9.29) 
6 

when ra^O. 

The case to = must be considered separately. It is 
formally obtained by setting 7 = in Eq. I|9.25|l . and 
Eq. H9.27|) shows that this amounts to letting a = 0. 
For m — 0, therefore, Eq. (|9.25(l reduces to Teukolsky's 
radial equation in Schwarzschild spacctime; the indepen- 
dent variable is now given by x = r/r + — 1. The rela- 
tion between ^> m (v ) 6) and the asymptotic value of tpo 
was already worked out, for Schwarzschild spacetime, in 
Sec. VIII E — Eq. P7) . There it was revealed that 
it is of the schematic form \f ,OT ~ r + ip (v,r ^> r + ,9,ip), 
and that it involves a derivative of the asymptotic field 
with respect to v. This relation is very different from 
what was anticipated in Eq. I|9.24[) . and therefore differ- 
ent from what is known to be true for to 7^ 0. These 
considerations imply that for to = and a ^ 0, the re- 
lation between *B m (v,0) and the asymptotic value of ipQ 
comes with an additional factor of M/1Z relative to terms 
with m/0. We conclude that it is appropriate to neglect 
the to = term in Eq. (|9.24|) . which becomes 

tpo = - A m \a m (v) + i(3 m (v)]x~ 2 (l + x)~ 2 

771^0 

x F(-4, 1; -1 + 2im 7 ; -x) 2 Y 2 m (6, ip), (9.30) 

where x = (r - r+)/( r + - r_), 7 = a/(r + - r_), A m 
is given by Eq. I|9.29|l . and a m (v), (3 m (v) are listed in 

Eqs. EHZMSSS- 

The functions * m (?j,0) are obtained by substituting 
Eq. (|9.3U|) into Eq. (|5.6|) and taking the limit r — ► r + , 
or x — * 0; we recall that %[)o(v,r,6,ip) is the Weyl 
scalar constructed with the Kinnersley tetrad, and that 
^(v, r + , 8, ip) is decomposed as in Eq. (|9.5() . Simple al- 
gebra, using a = %M, r± = M(l ± yT — x 2 ), and 
7 = |x(l -X 2 ) _1/2 , yields 

* m M) = - ^ ==2^ (1 + »m 7 )(l + 4m 2 7 2 ) 
12(1 + v/1 -x 2 ) 2 

x [a m ( W )+i/3 m («)] 2 y 2 m (^0). (9.31) 

This result holds when and it reveals that \l> m ^ 

is of order 1Z~ 2 ; as we have seen, when to = we have 
instead the Schwarzschild result *° = 0(M/TZ 3 ). 



D. Fluxes 

We now insert Eq. (|9.31l) into the approximate flux 
formulae of Eqs. l|9~TUj) and (|9~Tl) . The fact that *° = 
0(M/1Z 3 ) implies that the error term of Eq. I|9.10|l is in 
fact dominant, and we obtain 



This result indicates that to calculate (M) requires infor- 
mation that is not accessible to the leading-order analysis 
carried out here. To go beyond this leading-order calcu- 
lation should be feasible, but this lies beyond the scope 
of this work. 

A more definite result can be obtained for (J). Substi- 
tution of Eq. (|9.31() into Eq. I|9.11|l and integration over 9 
— recall the explicit forms of the spin-weighted spherical 
harmonics specified by Eqs. H9.14fl - 19.16fl — returns 



(j) 



M 5 X 
45 



£ [1 + (™ 2 



1) X 2 ] [4 + (to 2 - 4) X 2 



x (\a m (v) + i(3 m {v)\ 2 ) 



after simplification; notice that a factor of to 2 is canceled 
by the integral / \ 2 Y 2 m ^°(8, 0)| 2 sin (9 dO = 8/(5to 2 ). This 
becomes 



(j) = --M 5 x[(4-3x 2 )<£ 2 3+£f3 + Z? : 



#23) 



+ 4(l + 3x 2 )((£i 



4£ n 2 , 



22 ) 



after using Eqs. lj9~T7jl - l|9~2l|l . 

At this stage we introduce the invariants 



a0 



and 



El — £ a b£ a — £aji£ 

E 2 = £ ab s b £ a c s c = £ a ^£ a jS ^, 
E 3 = {£ ab s a s b ) 2 = {£ a ps a sP) 2 , 



Bt = B ab B ab = B af) B a ^ 

B 2 = B ab s b B a c s c =B apS PB a lS \ 

B 3 = (B ab s a S b ) 2 = (B a0 s a sP) 2 , 



(9.33) 
(9.34) 

(9.35) 

(9.36) 
(9.37) 

(9.38) 



where the unit vector s a gives the direction of the black 
hole's spin in the local asymptotic rest frame, and s a = 
s a e" is the corresponding spacetime vector. We therefore 
have J a = Js a , J = xM 2 , and (j) = (,j a )s a . In terms of 
these invariants we have, for example, £ 2 3 +£ 23 = E 2 — E 3 
and (£n - £ 22 ) 2 + 4£ 2 2 = 2E X - 4E 2 + E 3 . Our final 
expression for the rate of change of angular momentum 
is 



{j) 



:M 5 X 



(l + 3x 2 )(£i+Bi) 



45 

-3(4+17x 2 )(F 2 + J B 2 ) 
+ 15 X 2 (E 3 + B 3 



(9.39) 



(M) = 0(M 5 /K 5 



(9.32) 



This result reveals that (j) = 0(Af 5 /^ 4 )- 

The first law of black-hole mechanics implies that the 
rate of change of the horizon area is given by (k/8it)(A) = 
— ft-a(J) + 0(M 5 /H 5 ), with a leading term scaling as 
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M 4 /1Z 4 . The ratio Oh/k can be expressed in terms of M 
and \ = a/M, and we obtain 



(A) 



16?r M 5 X 2 



45 

~i{A + n x 2 ){E 2 + B 2 ) 

+ 15x 2 (E 3 +B 3 



(9.40) 



this scales as M 5 /7Z 4 . This result is equivalent to 
Eq. (3.5.39) of the book by D'Eath [l9|. 



E. Comparison with Thorne and Hartle 

The rate of change of angular momentum for a gen- 
eral body interacting with a tidal gravitational field was 
calculated by Thorne and Hartle |33 and their result dis- 
played in Eq. I|8.42[l . We wish to compare this general 
expression with our result for (J) displayed in Eq. 1)9. 39[) : 
recall that (J) = (J a )s a , with s a giving the direction of 
the angular-momentum vector. 

The quadrupole moments of a Kerr black hole im- 
mersed in an external universe include an intrinsic com- 
ponent that would be present even if the black hole were 
isolated, and an induced component that comes from the 
hole's tidal distortion. We write 



M ab = M a ' 



intrinsic 

t ^intrinsic _i 

'Jab — J a b 



M induced 
ab 



ab 



(9.41) 
(9.42) 



and we know that |3S 



M, 



intrinsic 
ab 



l 



M 6 x (Sab - 3s a s b ), 



rintrinsic n 

J ab — u - 



(9.43) 

We recall that M ab is the hole's mass quadrupole mo- 
ment, while J ab is its current quadrupole moment; both 
tensors are symmetric and tracefree. We wish to see if 
we can determine M^ duccd , and establish com- 

patibility between Eq. I|8.42[) and l|9.39|l . 

It is easy to show, by substituting Eq. (|9.43() into 
Eq. (|8.42|) . that the coupling between the intrinsic mo- 
ments and the tidal gravitational fields does not affect 
the magnitude of the angular-momentum vector; the 
only effect is to produce a precession of J a described 
by j a = e a bc n b p J c , where n<p = -M X £ a b s b is the pre- 
cessional angular velocity. We conclude that only the in- 
duced moments will contribute to (J), and we now seek 
to determine them. 

To ease the comparison between Eq. (I8.42|) and (|9.39() 
we set s° = (0,0,1) and compute (J) = (J 3 ), which 
we compare with the result displayed immediately be- 
fore Eq. !PH5|1 . This reveals that M^ duccd is partially 
determined by the relations M u -M 2 2 oc 16(1 + 3x 2 )£i2, 
M 12 oc -4(1 + 3 X 2 )(£ n - £22), M13 oc (4 - 3 X 2 )£ 23 , and 
M23 oc — (4 — 3x 2 )^23, where the (unique) constant of 
proportionality is equal to ^M 5 x- Analogous relations 



link § J^ ducod to B ab . These relations determine M™ duced 
(and J^ duced ) up to a term proportional to 5 a b — 3s a s b ; 
the coefficient must be a scalar formed from £ a b (or B a b), 
S a b, s a , and e abc , and the only possible candidate is an 
arbitrary function of x multiplying £ ab s a s b (or B ab s a s b ). 
We therefore arrive at 



M. 



induced 
ab 



45 



M 5 X 



Hx)(Sab - 3s a s b )£ cd s c s d 
+ 8(l + 3 X 2 )£ c {a e b)cd s d 
+ 30x 2 s (a e b)cd s c £ d e s e 

for the mass quadrupole moment, and 



(9.44) 



J. 



induced 



1 



-M 5 X 



c „d 



Kx){S a b - is a s b )B cd s c s 
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+ 8{l + 3 X 2 )B c {a e b)cd s d 
+ 30x 2 s {a e b)cd s c B d e s 



(9.45) 



for the current quadrupole moment, where A(x) and fJ>(x) 
are unknown functions of the hole's rotational parameter. 
We conclude that our results are indeed compatible with 
the general results of Thorne and Hartle • We observe 
that the relationships between the induced moments and 
the tidal gravitational fields have the schematic form 

M inducod _ M 5g t Jinducod _ M 5g ^ w j th & comp \i_ 

cated tensorial structure) , and that the moments scale as 
M 3 (M/1Z) 2 . These results follow expectation, and they 
are markedly different from those of Sec. VIII G; recall 
that for a nonrotating black hole the relationships involve 
an additional factor of M and a time derivative. 

We emphasize that the induced moments are only par- 
tially determined: the functions A(x) and n(x) cannot be 
determined by the comparison with Thorne and Hartle, 
because the terms to which they belong in M ab and J ab 
do not affect the magnitude of the angular-momentum 
vector. They produce instead a small fractional correc- 
tion of order (M/1Z) 2 to Qp, the precessional angular 
velocity. 



F. Comparison between Kerr and Schwarzschild 
results 



The main results of Sec. VIII, E qs. and jOty . 

or equivalently Eqs. (|8.40l) and l|8.41|l . hold to lead- 
ing order in M/TZ <C 1, and they reveal that for a 
Schwarzschild black hole, (M) = 0(M 6 /K 6 ) and ( j) = 
0(M 6 /1Z 5 ). On the other hand, the main results of this 
section, Eqs. (|9.32l) and l|9.39|l . hold to leading order in 
M/TZ <§; X: and they reveal that for a Kerr black hole, 
(M) = 0(M 5 /1Z 5 ) and (J) = 0(M 5 /1Z 4 ). The scalings 
are very different, and the condition M/TZ <§; x implies 
that the Schwarzschild results cannot straightforwardly 
be obtained from the Kerr results in a limit x — > 0. 

The origin of the difference in scalings can easily be 
understood in the special case of rigid rotation, for which 
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(M) and (j) are given by Eqs. l(P5|l - l|PT)) . 

(M) = fl(fl - Q K )K, (J) = {Q- Qh)IC, 

where K, is denned by Eq. (|4~32|) and O = 0(7£ _1 ) is the 
hole's angular velocity in the external spacetime. This 
argument was first presented to me by Kip Thorne (per- 
sonal communication), and it was then elaborated on by 
Alvi[3. 

Suppose, as we shall show below, that K, = 0(M 6 /K 4 ). 
In the case of nonrotating black hole we have Oh = 0, 
and it follows that (M) = Q 2 K, = 0{M 6 /1Z & ) and 
{j)=QX. = 0(M 6 /1Z 5 ); those are precisely the scal- 
ings obtained previously for a Schwarzschild black hole. 
The situation is different for a rotating black hole. In 
this case the condition M/TZ <C x implies that O <C Oh, 
and we have instead (M) = -00 H /C = 0(M 5 /1Z 5 ) and 
(J) = — S1h/C = 0(M 5 /1Z 4 ); those are precisely the scal- 
ings obtained previously for a Kerr black hole. Notice 
that in the case of Kerr, (M) and (J) are both pro- 
portional to Oh and therefore to x; this observation is 
confirmed by Eq. I|9.39(l . 

The different scalings reflect the different technical 
meanings assigned to the phrase "small- hole/slow- motion 
approximation:" For a Schwarzschild black hole we im- 
pose M/TZ <C 1 and we naturally have O ^ Oh; for a 
Kerr black hole we impose instead M/TZ <§; x and we con- 
sequently have O <C Oh- In generic situations (that is, in 
the absence of rigid rotation) the scaling argument given 
previously continues to apply, but ~ TZ^ 1 is now inter- 
preted as an inverse time scale associated with changes 
in the Weyl tensor of the external spacetime. 

An expression for K, can be obtained from the approx- 
imate relation (J) = — Oh/C which holds for a rotating 
black hole. We obtain 



hole is moving on a circular orbit, and we assume that 
the orbital angular momentum vector is either aligned or 
anti-aligned with the hole's spin vector: L ■ s = e = ±1. 
The hole's orbital angular velocity is then 







M + Mext 



b 3 



(9.47) 



where b is the orbital radius; the angular velocity is pos- 
itive when the orbital and spin angular momenta are 
aligned, and it is negative when they are anti-aligned. 

A simple calculation, along the lines of what was pre- 
sented in Sec. VIII H, yields 



8 i( M 



x(l + 3x*)V 1S (9.48) 



and 



M 



(J) ^ - r f ) (M + M cxtM l + 3x 2 )V- 

(9.49) 

where rj = MM cxt /(M + M cxt ) 2 is a dimensionless 
reduced-mass parameter and 



V = 



M + M e; 



« 1 



(9.50) 



is the relative orbital velocity. Notice that while we could 
not calculate (M) in the general case described in Sec. IX 
D, here it is simply given by (J) because the black hole 
is in rigid rotation around M oxt . The results of Eq. (|9.48|) 
and (|9.49|l agree with earlier expressions obtained by Alvi 
[18( ■ In the regime M ext ^ M they also agree with earlier 
results derived by Tagoshi, Mano, and Takasugi 



K = 



-M 6 (l + v/T^X 1 ) [8(1 + 3 X 2 ){E 1 + Bi) 



3(4 + 17x 2 )(£ 2 +£ 2 ) 



15x 2 (^3 + S 3 



(9.46) 



and we confirm that indeed, JC = 0(M & /1Z A ). Notice 
that there is no obstacle to taking the limit x ~~ * °f this 
expression. 



G. Black hole in a circular binary: Slow-motion 
approximation 

We now specialize the results of Sec. IX D to a Kerr 
black hole placed on a circular orbit in the weak gravi- 
tational field of an external body of mass M cxt . This is 
the slow-motion approximation, and we shall repeat here 
most of the steps described in Sec. VIII H. 

As before the tidal gravitational fields of the external 
universe are approximated by £ a b — &,ab and B a b — 0, 
where $ = — M CKt /\x — r cxt | is the Newtonian potential 
associated with the external body. As before the black 



H. Black hole in a circular binary: Small-hole 
approximation 

We now allow the Kerr black hole to move rapidly in 
the strong gravitational field of a Schwarzschild hole of 
mass M ext . We no longer restrict the size of V but we now 
impose M <C M oxt ; this is the small-hole approximation, 
and we shall repeat here most of the steps described in 
Sec. VIII I. 

Once more we take the orbit to be circular. In the 
standard Schwarzschild coordinates (t, r, 9, (f>) used in the 
background spacetime of the large black hole, the or- 
bital radius is b and the four-velocity of the small hole 
is u a = 7(1,0,0,0), where 7 = (1 - 3Af 0Xt /6)~ 1/2 is a 
normalization factor and 



I? 



(9.51) 



is the angular velocity; as in the preceding subsection 
e = ±1 gives the orientation of the orbital angular mo- 
mentum vector relative to the hole's spin vector, s a = 
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(0,0,-1/6,0). Calculation of £ a p, B a p using Eqs 
(|8.5|) . and substitution into Eqs. I|9.33f) - (|9.39[) gives 



(M) 



(i-2v 2 )(i-^v 2 ) 



(9.52) 



and 



(J) 



8/ M 



M extX (l + 3x 2 )^ 

O V Jw ext / 

(1-2^ 2 )(1-^^ 2 ) 
(1-3F 2 ) 2 



12 



(9.53) 



where V — \J Af oxt / b < 6 -1 / 2 is a measure of the hole's 
orbital velocity. These results agree with those of the pre- 
ceding subsection in a common domain of validity defined 
by M < M ex t and V < 1. To the best of my knowledge, 
the results of Eqs. (|9.47() and 19.48fl . complete with all- 
order relativistic corrections, have never appeared before 
in the literature. 



I. Black hole in a static tidal field 

For completeness we explore another special case of 
Eq. H9.39fl . in which the rotating black hole is at rest in a 
static tidal gravitational field. We wish to calculate the 
rate at which this black hole loses its angular momentum. 
This calculation was presented many times before, most 
notably by Hartle [5J, |55| , Tcukolskv [531 , Chrzanowski 
Thorne, Price, and Macdonald 12117 and Alvi (lgj . 
The point of this subsection is to illustrate how easily 
the classic spin-down result of Eq. (|9.55|) follows from 
Eq. 

We assume that the tidal gravitational field is purely 
electric in the local asymptotic rest frame of the black 
hole, and that it is axially symmetric in the arbitrary 
direction of the unit vector n a . With these specifications 
we have 

£ab = ~£(8 a b - 3n a n b ), B ab = 0, (9.54) 

where £ = £ ab n a n b . If, for example, the tidal field 
is produced by a body of mass Af cxt maintained at a 
fixed position r = bn relative to the black hole, then 
£ = — 2M 0Xt /6 3 . We assume that the black hole's an- 
gular momentum makes an angle a with respect to the 
direction of n a , so that s a n a — cos a. We then have 
£ a b sb = \£ (3cosan a — s a ), £ ab s a s b — i£(3cos 2 a — 1), 
and the invariants of Eqs. (|9.33|) - (|9.35|l are easily com- 
puted. After simplification we find that Eq. (|9.39|l re- 
duces to 



J 



;£ 2 M 5 X sm 2 a 



l-^(l-5sin 2 a)x 2 



, (9.55) 



which is the classic spin-down formula. 



Acknowledgments 

This work was supported by the Natural Sciences and 
Engineering Research Council of Canada. I am grateful 
to Eanna Flanagan, David Garfinkle, Karl Martel, and 
Clifford Will for helpful discussions. 



Note added 

Conversations with Kip Thorne and John Friedman 
(whom I thank) made me understand that the discussion 
of induced quadrupole moments inserted in Sees. VIII 
G and IX E is incomplete. I should have realized that 
the tidal-heating formulae of Eqs. 18.42|l and l|8.43|l al- 
low the determination of M ab only up to a term pro- 
portional to £ ab , and the determination of J ab up to 
a term proportional to B ab . Such terms do not par- 
ticipate in the tidal heating and leave (M) and (J a ) 
unchanged. It is therefore possible for a tidally dis- 
torted Schwarzschild black hole to have quadrupole mo- 
ments given by M ab = aM 5 £ ab + (32/45)M 6 £ ab and 
J ab = a'Al 5 B ab + {&/15)M 6 B ab , where a and a' are unde- 
termined dimcnsionless constants. Such moments would 
scale as M 5 /1Z 2 , which is the expected scaling. 



APPENDIX: PERTURBATIONS OF A 
SCHWARZSCHILD BLACK HOLE 

In this Appendix I collect a few key results from the 
theory of gravitational perturbations of a Schwarzschild 
black hole 0, 0, |H, 0] . I employ a covariant/gauge- 
invariant formalism that was inspired by the work of 
Gerlach & Sengupta [HIH and Sarbach & Tiglio [H. 
These results are presented without derivation; details 
can be found in Martel's PhD dissertation [Toj . 



1. Background metric 

The Schwarzschild metric is expressed as 

ds 2 = g tj dx l dx j + r 2 n AB d6 A dO B , (A.l) 

in a form that is covariant under two-dimensional coordi- 
nate transformations x l — > x l . The indices i, j, fc, . . . run 
over the values and 1, and the indices A, B, C, . . . run 
over the values 2 and 3. The traditional Schwarzschild 
coordinates are x l = (t, r), and in the text we use the 
ingoing Eddington-Finkclstein coordinates x 1 — (v,r), 
where v = t + r + 2Mln(r/2M - 1), with M den oting 
the black- hole mass. In the metric of Eq. (|A.1(1 . r is 
viewed as a scalar function of the arbitrary coordinates 
x l , and Qab = diag(l, sin 2 9) is the metric on the unit 
two-sphere. 

We use gij and its inverse to lower and raise all lower- 
case Latin indices. And in this Appendix, contrary to 
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previous usage in the body of the paper, we use Qab 
and its inverse to lower and raise all upper-case Latin 
indices. We indicate covariant differentiation with re- 
spect to a connection compatible with gij with a dot: 
9ij.k = 0. And we indicate covariant differentiation with 
respect to a connection compatible with £Iab with a 
colon: VIab-.c = 0. 



2. Spherical harmonics 

The tcnsorial nature of the spherical harmonics refers 
to the unit two-sphere, whose metric is VLab- The def- 
initions adopted below agree with those of Regge and 
Wheeler The Levi-Civita tensor on the unit two- 
sphere is denoted sab, with ee^ = siylO. 

The scalar harmonics are the usual spherical-harmonic 
functions Y lm (9 ), which satisfy the eigenvalue equation 
tt AB Y% + 1(1 + l)Y lm = 0. 

Vectorial spherical harmonics come in two types. The 
even-parity harmonics are 



and 



Y : ™ (even parity), 



while the odd-parity harmonics are 

(odd parity). 



vim 

X A - 



,_ B^rlm 
~ E A 1 :B 



(A.2) 



(A.3) 



The vectorial harmonics satisfy the following orthogonal- 
ity relations: 



y ( a y^' dn = I xtx l A m ' dn = i(i + 1) s w s r} 



and 



Y^x i ; m ' dn = o, 



(A.4) 



(A.5) 



where an overbar indicates complex conjugation and 
dft = sin 8 d8d(f>. 

Tensorial spherical harmonics come in the same two 
types. The even-parity harmonics are 



Y lm n AB , 



Y 



:AB 



(even parity). 



It is useful to define also 



A AB — 



Y l % B + ±i(i + i)Y lm n AB . 



(A.6) 



(A.7) 



By virtue of the eigenvalue equation for the scalar har- 
monics, tt AB Z 1 ab = 0; these harmonics are therefore 
tracefree. The odd-parity harmonics are 



Xlra 
AB 



viva 
~^{A:B) 



(odd parity); 



(A. 



these are also tracefree: il X 



AB ylm 
AB 



0. We record the 



following relations between the tensorial harmonics and 
the spherical harmonics of spin- weight s = ±2 |60|| : 



^AB 



1 



x (_ 2 Y lm e A e B + 2Y lm e A e B ) (A.9) 



vim 

A iB — 



-V(Z-1W + !)(* + 2) 
x(_ 2 Y lm e A e B - 2 Y lm e A e B ), (A.fO) 



where the vectors ea = (1, i sin9)/^/2 satisfy Q, AB eA^B = 
n AB cacb = and U, AB e A ^B = 1- The tensorial harmon- 
ics satisfy the following orthogonality relations: 



x AB x l ^'dn 



-a-i)i(i + i)(i + 2)s w s ri 



and 



~,AB V V 



x A g dn = o. 



(A.ll) 



(A.I2) 



The spherical-harmonic functions constructed thus far 
are complex, and they are all proportional to e™^. We 
shall also need a set of real spherical harmonics of degree 
1 = 2, which we denote Y m , Y™, XJ, X™ B , ZJ B , and 
X™ B . The label m runs over the set {0, lc, Is, 2c, 2s}; the 
numerical part of this label refers to the azimuthal index 
to, and the letter indicates whether the corresponding 
scalar-harmonic function is proportional to cos(to^) or 
sin(TO0). Explicitly, 



Y° 




cos 2 6>- 1), 


Y lc 


= sin 6 


'cos#cos</>, 


Y ls 


— sin 6 


1 cos 9 sin cj), 


Y 2c 


= sin 2 


8cos2cj), 


y2s 


. 9 
= Sill 


6>sin2</>. 



(A.13) 



The vectorial and tensorial harmonics are generated by 
acting on Y m with the same differential operators as 
those involved in Eqs. I|A.2|) . I|A.3|) . (fA.Bfl — (fA.Sfl . 



3. Odd-parity perturbations 

The odd-parity perturbations of the Schwarzschild 
metric are those which are expanded in terms of odd- 
parity spherical harmonics. This sector of the metric 
perturbation is given by 

6g iA (x\9 A ) = h l (x l )X l X l (8 A ), (A.14) 
5g AB (x\9 A ) = foix^Xfy^). (A.15) 

We suppress usage of the Im label on the fields hi and 
/i2, and it is understood that the right-hand sides are 
summed over I and to. It can be shown that the combi- 
nations 



1 



r.iho 



(A.16) 
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are invariant under odd-parity gauge transformations. 
The linearized Einstein field equations are then natu- 
rally expressed in terms of hi and its covariant deriva- 
tives. One of these equations is required in the text: In 
the absence of sources it can be shown that 



0. 



(A.17) 



The remaining field equations can be manipulated to 
form a one-dimensional wave equation for the master 
variable 

*rw = -r'%, (A.18) 
r 

which is evidently gauge invariant. The function 
$rw(^') is known as the Regge- Wheeler function Q , and 
in the absence of sources it satisfies the differential equa- 
tion 

'1(1 + 1) 6M~ 



RW 



0. 



(A.19) 



where □'J = g^^.ij is the one-dimensional wave op- 
erator acting on the scalar function ^>(x l ). It is well 
understood that in a specified gauge, all components of 
the odd-parity metric perturbation can be reconstructed 
from the Regge- Wheeler function. 



4. Even-parity perturbations 

The even-parity perturbations are expanded in terms 
of even-parity spherical harmonics. This sector of the 
metric perturbation is given by 



Sg i:j (x l 
&giA{x l 
SgAB(x l 



3l (x l )Y^(9 A ), 



(A.20) 
(A.21) 



l ) = r 2 [K{x l )Y lm (0 A )n 



AB 



+ G(x l )Y%(9 A )]. (A.22) 

Once more we suppress usage of the Im label on the fields 
hij, ji, K, and G, and it is understood that the right- 
hand sides are summed over I and m. The combinations 



hij — hij 



2e 



(i-3) 



K = K r 

r 



(A.23) 



where = ji — ^r 2 G_i, are invariant under even-parity 
gauge transformations. The linearized Einstein field 
equations are then naturally expressed in terms of these 
fields and their covariant derivatives. They can be ma- 
nipulated to form a one-dimensional wave equation for 
the master variable 



2r 



1(1 



1) 

where A = (/ — l)(l + 2) + 6M/r. The function ^zm(x 1 ) 
is known as the Zerilli-Moncricf function [J, |43| , and it is 
evidently gauge-invariant; it satisfies a differon! ial oqua- 



K 



j(r^hij 



(A.24) 



tion similar to Eq. l|A.19fl . but with a more complicated 
potential. The normalization of the Zerilli-Moncrief func- 
tion is chosen so as t o ag ree with the definition proposed 
by Lousto and Price |13| . 



5. Waveforms and energy radiated at infinity 

When examined near future null infinity, the gravi- 
tational perturbations of Eqs. (|A.14f) . (|A.15fl . (jA.20|) - 

(|A.22|) can be presented in an outgoing-radiation gauge 
that permits an easy identification of the radiative field. 
It can be shown that the two fundamental polarizations 
of the gravitational waves are given by 



h+-ih x = iy;v(i-iw + i)(i+2) 

It * — ' 



l in 

2i 



^(u')du' 



(A.25) 



2M\n(r/2M - 1) is retarded time, 



\v A ) are spherical harmonics of spin- weight 
The fact that the waveforms are expressed 



where u = t 
and _ 2 Y lm (6 
s = -2 

in terms of an integral of the Regge- Wheeler function 
means that this master variable is rather ill-suited to de- 
scribe the radiative aspects of the metric perturbation. 
An alternative choice of master variable, which is free of 
this blemish, was proposed by Cunningham, Price, and 
Moncrief |6ll |62| ; it was recently revived by Jhingan and 
Tanaka [63^ 

The energy and angular momentum radiated to infinity 
are given by 



(E) 



(J) 



g L £(i +!)(/ + 2) 

Ira 

x(4|^( U )| 2 + |^™(u)| 2 ), (A.26) 

-L^-l)Z(i+l)(Z + 2)(im) 
Im 

x (4*j&(«) f U ¥^(u')du' 

(A.27) 



*'zm(«)*£m(«) 



The expressions are very similar to the horizon-flux for- 
mulae of Eqs. (jTSj) and (J73- 
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